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The  motion  of  a  thin  wall  is  treated  in  general  relativity  in  the 
case  where  the  regions  on  either  side  of  the  wall  have  non-zero 
stress-energy.  The  space-times  studied  here  are  the  following:  (i) 
spherically  symmetric  with  vacuum  energy,  (ii)  plane-symmetric  with 
vacuum  energy,  (iii)  plane-symmetric  fluids  with  pressure  equal  to 
energy  density,  and  (iv)  plane-symmetric  dust  with  vacuum  energy. 
In  connection  with  the  latter  case,  new  classes  of  exact  solutions  to 
Einstein's  equation  are  presented,  corresponding  to  plane-symmetric 
non-homogeneous  dust  universes.  The  equation  of  state  of  the  thin 
wall  is  taken  to  be  x=To,  where    x      is  the    tension    in  the  wall,    o    is  the 


energy  density,  and  0<_  r<_l.  Particular  cases  of  interest  are  r  =  0, 
corresponding  to  a  dust  wall,  and  r=l,  corresponding  to  a  domain  wall. 
The  dynamics  of  the  thin  wall  depend  on  both  the  stress-energy  of 
the  wall  and  that  of  the  space-time  on  either  side  of  the  wall. 
Oscillation  about  a  point  of  stable  equilibrium  is  found  to  be  possible 
for  walls  with  certain  values  of  the  parameters  in  cases  (i)-(iii).  In 
each  such  case,  the  metric  on  either  side  of  the  wall  is  static.  There 
are  various  other  possible  dynamics,  such  as  expansion  followed  by 
collapse,  deflation  followed  by  expansion,  simple  translation,  and 
asymptotic  approach  to  points  of  unstable  equilibrium.  Finally,  the 
global  structure  for  each  type  of  space-time  is  discussed.  This  permits 
the  elucidation   of  the   structure   of  the   corresponding   wall   space-times. 


vn 


CHAPTER  1 
PRELIMINARIES 

Introduction 


One  of  the  intriguing  possible  consequences  of  the  hot  big  bang 
model  of  cosmology  is  that  the  universe  may  have  undergone  phase 
transitions   in   its  early   history.      These   phase   transitions   can   give  rise 

to   various   soliton-like   structures,   such   as   magnetic   monopoles,   cosmic 

1  ? 
strings,    and   domain   walls.1 '^   The  dynamics   of  the   latter  hypothetical 

object,  as  a  special  case  of  a  more  general  class  of  thin  walls,  will  be 
the  focus  here.  The  space-time  stress-energy  on  either  side  of  the 
wall  will  be  taken  to  contain,  in  turn,  (i)  vacuum  energy,  (ii)  fluid  with 
pressure  equal  to  energy  density,  and  (iii)  dust  with  vacuum  energy. 
In  connection  with  the  last  case,  new  classes  of  exact  solutions  to 
Einstein's  equation  will  be  presented,  corresponding  to 
non-homogeneous  dust  universes.  Most  previous  work,  with  the 
exceptions  noted  below,  has  dealt     simply  with  vacuum. 

Thin  walls  may  arise  from  contributions  by  particle  fields  to  the 
space-time  stress-energy.  A  simple  field  theory  giving  rise  to  thin 
walls  contains  a  real  scalar  field  §  with  a  potential  V(<|>),  where  V(<|>)  is 
taken  to  have  two  local  minima.  Let  <)> ,    and  <]>_  denote  the  values  of    § 

corresponding  to  these  minima.     During  a  phase  transition,     the  field     § 

1 


takes  on  values  near  ty+  in  one  region  of  space  (  the  "  §  +  -  domain"), 
and  near  §     in  an  adjacent  region  of  space  (  the  "  §_  -  domain   ").   The 

region  between  the  §+  and  $     domains  has  a  high  energy  density,  both 

because  the  field  takes  on  values  which  don't  minimize  V,  and 
because  there  is  energy  associated  with  spatial  changes  in  the  field. 
This   high-energy  region  is  designated   a  domain  wall. 

To    make    these    ideas    more    concrete,    take    a    particle    field    § 

together       with       a       potential       V(())),       and       the       Lagrangian 

9 
L=-(l/2)(Va<|))    -V((j)).    The    contribution    to    the    energy-momentum    of 

space-time  by  this  field  is  given  by 

Tab  =  Ca^b*)  "  (l/2)gab(vC(t)Vc^+  2Vfo»-  0-1) 

If  <\>  is  approximately  constant  throughout  a  region  of  space-time,  then 
the  first  two  terms  on  the  right  of  (1.1)  are  negligible  in  that  region, 
whereas    the    third    term    is    a    multiple    of    the    metric.    Hence    in    the 

4>+ -domain,   say,   the  contribution   to  the  stress-energy   is 

Tab+  =  -V(<Ugab    .  d-2) 

A  solution  to  Einstein's  equation  with  a  contribution  such  as  this  to 
the    stress    energy      is    the   same   as    a   solution    of  equations    having   a 


cosmological    constant,    with       A+  =  kV(<)>  +  )  ,  where  k    is    Einstein's 
constant.  An  identical  situation  holds  in  the  <J>  -domain. 

Under  the  conditions  described  above,  the  thin  wall 
approximation  developed  by  Israel  provides  a  possible  model  for  a 
space-time  containing  such  a  wall.  To  make  the  calculations  tractable, 
the  wall  is  assumed  to  have  plane  or  spherical  symmetry,  at  least, 
and  often  complete  reflection  symmetry  as  well.  Mathematically,  this 
latter    condition    is    expressed    by    Ka^  +  =  -Ka^_,    where    Ka^   is  the 

extrinsic   curvature   (defined   in   the   next   section).      The   work   done   on 
planar   walls   by   Vilenkin   ,  Ipser  and  Sikivie   ,  and     Ipser"    treated  the 

case   V((|)+)=V(<))  )=0.  Most  of  the  work  on  spherical  domain  walls  has 
been  restricted   to  special  values   of  the  four  parameters   A+  ,  A    ,  m+  , 

o    c    n     1  n  10 

and     m_   '    ''"      ,   with    some   numerical   work   on    the   general   case      . 

Garfinkle   and   Vuille  9   studied  the  general  case  V(<|> ,  )*V((|>_),  for  both 

spherical   and  plane-symmetric  walls.      Very  little  has  been  done  in  the 

case    where    there    is    energy-momentum    other    than    vacuum    energy. 

20 
Kuchar  w  examined  special  cases  involving  collapse  of  a  charged  shell, 

9  1 
while     Maeda         specialized    the    method    of    Israel    to    space-times 

containing   an   ambient  fluid   without  examining   wall   dynamics. 

To  find  the  motion   and  the  gravitational   field   of  the  wall,  it  is 

first    necessary    to    solve    Einstein's    equation    in    the    presence    of    the 

stress-energy    of    the    wall,    the    vacuum    energy    of    the    domains    on 


either  side  of  the  wall,  and  the  energy-momentum  contributed  by  any 
other  distributions  of  matter  or  fields.  This  gives  the  gravity  field  of 
the  space-time.  Israel's  thin  shell  approximation,  discussed  in  the  next 
section,  gives  an  equation  for  the  jump  in  the  extrinsic  curvature  in 
terms  of  the  stress-energy  of  the  fluid  composing  the  wall.  This  leads 
to  the  equation  of  motion  giving  the  dynamics  of  the  wall  in  the 
space-time.  In  Chapter  Two,  the  dynamics  of  thin  walls  in  regions 
with  vacuum  energy  will  be  addressed,  for  both  spherical  and  planar 
symmetry.  In  Chapters  Three  and  Four  the  space-times  under  study 
will  be  plane-symmetric  and  contain  fluid  with  pressure  equal  to 
energy   density,   and   dust   with   vacuum   energy,   respectively. 

2.  The  Thin  Shell  Approximation  of  Israel 

In  the  thin  shell  approximation  of  Israel  a  thin  shell  of  matter 
is  idealized  as  having  zero  thickness  (for  a  rigorous  justification  of  this 
approximation,  see  reference  22).  The  history  of  the  thin  shell  is  just  a 
three-dimensional  timelike  hypersurface,  and  thus  the  space-time 
consists    of    a    smooth    manifold    M,    a    metric    ga^    ,    and    a    smooth 

three-dimensional  timelike  hypersurface  Z.  The  hypersurface  divides 
the  manifold  M  into  two  manifolds  with  boundary,  called  M+  and  M". 
In    M     and    M"    separately    g   l    is    a    smooth    solution    of    Einstein's 

equation.  Let  na  be  the  unit  normal  to  Z  pointing  from  M"  to  M+.  Then 


in  each  region   M+   and  M"   the   space-time   metric   gab   induces  on   the 

hypersurface    2   an  intrinsic   metric   h  b   and  an  extrinsic  curvature  K   u 
given    by 


hab=    §ab  -nanb  O-3) 

:ab  =  haCvcnb 


Kah=hflcVrnh     .  (1.4) 


The    latter    quantity    may    differ    depending    on    whether    it    is 
calculated    in    M+  or   M".   Denote  by  6+   the  value  of  a  tensor  0  on  Z 

calculated    in    M+;  define  G_  similarly.    The    appropriate    conditions    on 

the  change  in  hab  and  Kab  across    I    are3 

hab+-hab-=°  d-5) 

Kab+  -  Kab-  =  -*[Sab-(l/2)  S  hab]  (1.6) 

where    Sab    is    the   so-called    surface    stress-energy    tensor   of   the    shell 
and  the  scalar  S  is  given  by  S=habSab.    Sab     may  be  interpreted  as  the 

limit,    as    the    thickness    of    the    shell    goes    to    zero,    of    the    shell 
stress-energy   integrated   over   the   thickness    of  the   shell. 

Equations    (1.5)    and    (1.6),    together    with    Einstein's    equation    in 
M       and  M"    separately,   are  equivalent  to  Einstein's  equation   in   all   of 


M.  Thus  one  can  obtain  a  solution  to  Einstein's  equation  in  the  thin 
shell  approximation  as  follows:  first  choose  two  smooth  solutions  to 
Einstein's  equation  on  the  two  manifolds  M+  and  M".  Then  in  each 
manifold  choose  a  timelike  hypersurface  (Z+  and  £"  respectively). 
The  hypersurface  Z+  divides  M_+  into  two  regions.  Choose  one  of  the 
two  regions.  Then  the  union  of  this  region  with  Z+  forms  a  manifold 
with  boundary  which  we  call  M+.  In  the  corresponding  way  obtain  the 
manifold  with  boundary  M"  from  IT  and  one  of  the  two  regions  of 
M ".  Then  identify  the  points  of  Z+  and  Z"  to  produce  a  single 
hypersurface  £  ,  and  thus  a  single  manifold  M  which  is  essentially  M+ 
and  M"  "glued  together."  The  surfaces  Z+  and  S"  and  the  map  which 
identifies  them  must  be  chosen  so  that  equations  (1.5)  and  (1.6)  are 
satisfied.  When  the  surfaces  and  the  map  are  so  chosen,  the  resulting 
space-time  M  is  a  solution  of  Einstein's  equation  in  the  thin  shell 
approximation. 

The     stress-energy     tensor     of    the     general     (two-dimensional) 
perfect  fluid  is  given  by 

Sab  =(a-T)vavb  -xhab  (1/7) 

where  a  and  x  are  scalars,  the  energy  density  and  tension  of  the  wall, 
respectively,  and  va  is  a  unit  timelike  vector,  the  space-time  velocity 
of  the  fluid.  The  properties  of  a  perfect  fluid  are  given  by  its  equation 


of  state,  a  relation  which  gives  x  as  a  function  of  o  .  A  domain  wall 
has  a  particularly  simple  equation  of  state:  x  =  a.  It  then  follows  from 
conservation   of  stress-energy   that      o    is  constant     and  that 

Sab=-°hab      ■  d-8) 

In  this  work  a  more  general  equation  of  state,  x  =  To,  will  be 
assumed,  where  T  is  a  constant  satisfying  0  <  r  <  1  .  This  equation  of 
state  is  also  used  in  references  5,6,7,  and  18.  Here,  only  the  case  a  >  0 
is  considered.  Call  the  class  of  thin  perfect  fluid  walls  with  this 
equation  of  state  T -walls".  Domain  walls  correspond  to  T-walls  with 
r=l.  Hence  the  class  of  T-walls  is  more  general  than  domain  walls,  but 
less  general  than  the  whole  class  of  perfect  fluid  walls  with  arbitrary 
equation  of  state.  It  follows  from  the  equation  of  state  that  the 
stress-energy  of  a  T-wall   is 

sab  =  o[(i-r)vavb-rhab]   .  (i.9) 

Using  equation  (1.9)  in  (1.6)  gives  the  change,  for  a  T-wall,  in  the 
extrinsic   curvature  K&b  across  S  : 

Kab+  -  Kab-  =  -Ka{(l-Dvavb  +  (l/2)habJ     .  (1.10) 


This    is    the    key    equation    that    will    be    used    to    obtain    equations    of 
motion  for  thin  walls  under  the  various  conditions.  As  will  be  seen  in 
the  next  section,   the  left-hand   side  involves   components   of  the   metric 
and    its   derivatives,    and    hence    depends    on    the    ambient    stress-energy 
on  either  side  of  the  wall. 


CHAPTER  2 
THIN  WALLS  IN  REGIONS  WITH  VACUUM  ENERGY 


Spherical  Walls 

To  find  the  motion  and  gravitational  field  of  a  spherical  wall  it  is 
first  necessary  to  find  two  spherically  symmetric  solutions  of  the 
vacuum  Einstein  equation  with  cosmological  constant.  A  space-time  is 
spherically  symmetric  if  its  symmetry  group  contains  the  symmetry 
group  of  the  2-sphere.  The  general  spherically  symmetric  solution  of 
the  vacuum  Einstein  equation  with  cosmological  constant  A  is  either 
the  Schwarzschild-de  Sitter  metric  or  the  Nariai  metric.23  In  the 
Nariai  space-time  all  the  2-spheres  of  symmetry  have  the  same  area, 
and  so  it  is  unsuitable  for  describing  an  expanding  or  contracting  shell 
and  will  not  be  considered  here.  The  Schwarzschild-de  Sitter  metric 
has  the  following  form.    3 

ds2  =  -Fdt2  +  F_1dr2  +  r^de2  +  sin20  d(J>2)  ,  (2.1) 

where     the  function  F  is  defined  by 


F=    1  -  2mr_1     -  (l/3)Ar2  ,  (2.2) 


10 

and  m  is  a  constant.  This  space-time  can  be  regarded  as  a  black  hole 
of  mass  m  in  a  space-time  which  is  asymptotically  de  Sitter.  In  the 
case  m=0  the  space-time  is  de  Sitter  space-time,  while  in  the  case  A=0 
the  space-time  is  Schwarzschild.  Note  that  the  expression  in  equation 
(2.1)  is  badly  behaved  at  points  where  F=0.  This  is  simply  because  the 
coordinate  t  is  badly  behaved  at  these  points.  In  fact,  the  metric  is 
smooth  for  all  values  of    r     greater  than  zero. 

Each  of  the  space-times  M+  and  M_~  are  taken  to  have  a  metric  of 
the  form  given  in  equation  (2.1).  Let  A+  and  m+  be  the  values  of  the 
constants  A  and  m,  respectively,  for  the  metric  on  M  .  and 
correspondingly    define     A     and  m  .    The    cases    considered    here    are 

those  for  which  all  four  of  these  parameters  are  greater  than  or  equal 
to  zero  . 

Next,  it  is  necessary  to  find  timelike  hypersurfaces,  X+  and  £",  in 
M  and  M".  and  calculate  the  intrinsic  metric  and  extrinsic  curvature 
induced  on  these  hypersurfaces.  Since  the  wall  is  spherically 
symmetric,  £  and  S"  must  be  spherically  symmetric  also.  First 
consider  a  general  spherically  symmetric  timelike  hypersurface  I  in 
a  space-time  M_  of  the  form  given  in  (2.1).  Later,  these  results  will 
be  applied  to  M_+  and  M  "  separately.  A  three-dimensional 
hypersurface  in  a  four-dimensional  space-time  is  specified  by  giving 
the  four  space-time  coordinates  as  functions  of  any  three  parameters. 
Choose   the  coordinates   r,   0,  and  §   as   the  three  parameters.   Then   the 


1 1 

surface  is   specified  by  giving  the  remaining  coordinate  t  as  a  function 

of    the    three    parameters.        Since    the     surface   £_  is  spherically 

symmetric,     t  on     I_     must  be  a  function  of     r     alone.  Thus     I    is 
specified  by  a  relation  of  the  form 

t  =  T(r)  .  (2.3) 

A   normal  co-vector  N&  to  the  hypersurface     X    is  given  by 

Na  =  Va(t  -  T(r))    =    (  1  ,  -dT/dr  ,0,0).         (2.4) 

Define  the  scalar  Q  by  Q=NaN     .  Then  equations  (2.1)  and  (2.4)  give 


Q  =  F(dT/dr)2  -F"1  .  (2.5) 


Since  Z  is  timelike,  any  normal  to  I  must  be  spacelike.  Thus  Q  must 
be  greater  than  zero.  Thus  the  general  timelike  spherically  symmetric 
hypersurface  in  M.  is  given  by  a  relation  of  the  form  in  equation  (2.3) 
satisfying  the  condition  that  the  right-hand  side  of  (2.5)  is  greater 
than  zero.  Let  na  be  a  unit  normal  to  I.  Then  it  follows  from  the 
definition  of  Q  and  Na  that 
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na  =  eQ-1/2Na,  (2.6) 


where     e  =  ±l.     Note     that     there     are     two     different     unit     normals 
corresponding  to  the  two  different  possible  values  of  e. 
Define  the  tensor  q  l  by 

qabdxadxb  =  r2(d92  +  sin20  d<))2)    .  (2.7) 

Then  the  metric     can  be  written  as 

ds2  =  -Fdt2  +  F_1dr2  +  qabdxadxb  .  (2.8) 

Thus  qab  is  the  "2-sphere  part"  of  gab.  A  perfect  fluid  wall  has  a  fluid 

velocity  va  which  is  a  timelike  unit  vector,  tangent  to  the  wall.  A 
spherically  symmetric  wall  must  have  spherically  symmetric  va.  Thus 
the  conditions  on  va  are 

vava  =  -1  ,  (2.9a) 

vana  =  0  ,  (2.9b) 

vaqab  =  0  .  (2.9c) 

Using  these  conditions  and  equations  (2.4),  (2.6),  and  (2.7),  obtain 


13 

va  =  5Q"1/2(  -F(dT/dr),     F'1  ,  0  ,  0  )  ,  (2.10) 

where    8=±1.  It     follows,  using  equations  (2.4)-(2.8)     and  (2.10),  that 

Sab  =nanb  "vavb  +  %b  ■  (211> 

Inserting  equation   (2.11)  into   (1.3)  yields   the   following   expression   for 
the   intrinsic    metric: 

hab  =  "vavb  +%b  ■  (2-12) 

Since   hab    is    spherically    symmetric,    it    determines    a    set    of 

2-sphere   subspaces   of  L  (the  spaces  of  constant  r).  The  metric  induced 

on  a  constant  r  subspace  is  qab  and  the  area  of  the  subspace  is  4kt   . 

Thus  both     qak     and    r    are  determined  by     h  u    .  It  then  follows  from 

equation  (2.12)  that     v&vj)     is  determined  by    h^. 

Using    equations    (1.3)    and    (1.4)    and    the    fact    that    na    is 
hypersurface   orthogonal,   it  can   be   shown   that 

Kab  =  (l/2)£nhab  ,  (2.13) 
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where     £n     denotes     the  Lie  derivative  with  respect  to  the  vector     nJ 
Inserting   equation    (2.12)   results   in 

Kab   =  (l/2)(£nqab  -  va£nvb  -  vb£nva)  .  (2.14) 


The  Lie  derivatives  of     qab    and    va      are  computed   by  expressing  the 
Lie    derivative    in    terms    of    the    coordinate    derivative    operator24  3 
and  using  the  expressions  for  the  coordinate  components  of     q  b  ,  v    , 
and    na.    Some   straightforward   algebra   then   yields 

Mab  =  ^c^ab  +   %bd2Ln°  +  %cabnC= 

=  -2er-1Q-1/2F(dT/dr)qab  ,  (2.15) 

£nva  =  nbabva  +  Va"b= 

=  -(d/dr){eQ-1/2F(dT/dr)}va  .  (2.16) 

Define  the  scalar  L  by 

L=  eQ"1/2F(dT/dr)    .  (2.17) 

Then      equations   (2.14)-(2.17)   give 
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Kab  =  (dL/dr^Vb-r^Lq^    .  (2.18) 

The    quantity    L    can    be    expressed    in    two    different    and    useful 
ways.  First,  using  equations  (2.4),  (2.6).  (2.8)  and  (2.17),     obtain 

naVar  =  eQ"1/2NaVar  =  -L  =  -sILI  ,  (2.19) 

with  the  quantity     s     defined  by 

s  =  L/ILI   .  (2.20) 

The  quantity  s  can  be  interpreted  as  follows:  the  surface  E  divides 
M  into  two  regions.  A  region  is  an  "exterior"  if  r  increases  along  the 
direction  normal  to  E  and  pointing  into  the  region.  A  region  is  called 
an  "interior"  if  r  decreases  along  the  direction  normal  to  £_  and 
pointing  into  the  region.  Then  L  divides  M_  into  an  exterior  and  an 
interior.  If  s  =  1  then  na  points  to  the  interior.  If  s  =  -1,  then  na 
points   into   the   exterior. 

The     other     useful     expression     for     L     can     be     derived     using 
equations  (2.1),  (2.11),  and  (2.19)   : 

F  =  gn  =  (n1)2-(v1)2  =  L2  -  (v1)2  .  (2.21) 
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Then   equation   (2.20)  results   in 

L^sfF  +  Cv1)2}1/2.  (2.22) 

Choose  a  constant  r  subspace  S  of  Z  and  for  each  point  p  e  X 
and  let  X(p)  be  the  proper  time  along  the  integral  curve  of  va  from 
S    to    p.  Then  there  is  a  function    R(k)    such  that  on    £ 

r  =  R(X).  (2.23) 

It  follows  that  on    Z 


vavar=  v1   =  dR/cR  .  (2.24) 


Thus   from   equation   (2.22)   there  results 

L  =  s(F  +  (dRAR)2)1/2.  (2.25) 

Next  it  will  be  shown  that  the  surface  E_  is  essentially 
determined  by  the  function  R(k).  Given  the  function  R(X.)  and 
equation  (2.23),  the  quantities  dR/dX  and  L  are  determined  as 
functions  of  r.  Equations  (2.3)     and  (2.24)  give 
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vaVat   =   (dT/dr)(dR/cR)  .  (2.26) 


Then     equations  (2.8),  (2.9a),  (2.24),  and  (2.26)  result  in 

-1  =  gabvavb  =  -F(vaVat)2  +  F^v^r)2  = 

=  (dR/d?i)2{-F(dT/dr)2  +  F1 }  .  (2.27) 

Solving    for   (dT/dr)^   and   using  equation   (3.25)   yields 

(dT/dr)2  =   (L/F(dR/d?i)}2  .  (2.28) 

Since  the  right-hand  side  of  (2.28)  is  determined  as  a  function  of  r,  it 
follows  that  the  function  T(r)  is  determined  up  to  a  change  of  sign  and 
addition  of  a  constant.  Thus  if  Z.  and  Z.'  are  two  surfaces  in  M 
corresponding  to  the  same  R(^.),  then  one  can  take  Z_  to  Z.'  by 
changing  the  coordinate  t  on  M.  by  the  addition  of  a  constant  and 
(possibly)  a  change  of  sign.  But  this  coordinate  transformation  is  an 
isometry.  Thus  the  space-time  M.  with  the  surface  Z.  is  identical  to 
the  space-time  M_  with  the  surface  Z'.  So  in  this  sense  the  surface  Z 
is  determined  by  the  function  R(X). 

These   results   can   now   be   applied   to   the    surfaces    Z+    and    Z", 
which  are  glued  together  to  form  the  surface     Z.  Since    ha^+  -  ha^_=  0, 
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it  follows  that  if  0  is  any  tensor  defined  on  Z+  and  Z  such  that  0 
is  determined  by     h^,  then    0+  -  0    =  0.  Thus 

W  -  lab-  =  °  >  (2-29) 

r+  -  r_  =  0  ,  (2.30) 

(vavb>+-  (vaV-  =  °  •  (2-31) 

Choose  the  origin  of  X,  and  the  direction  of  va  to  be  the  same  on  Z+ 
and    I".  Then  it  follows  from  (2.31)  that 

X+-X_  =0,  (2.32) 

and  from  equations   (2.23),   (2.30),  and   (2.32)  that 

R(X+)  -  R(X_)  =  0  ,  (2.33) 

so  that  the  surfaces  Z+  and  £"  are  described  by  the  same  function 
R(X).  So  when  the  points  of  Z+  are  identified  with  those  of  Z",  each 
point  must  be  identified  with  a  point  which  has  the  same  value  of  X  in 
order  to  comply  with  equation  (2.30).  Similarly  each  point  must  be 
identified  with  a  point  which  has  the  same  values  of  6  and  $  (up  to  an 
overall    rotation)   in   order   to   comply   with    equation    (2.29).    When    the 
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surfaces   are  identified   in   this   manner,   it  follows  from  equation   (2.12) 
that  the  condition     hab+-hab_=0     *s  satlsfied- 

Now  the  only  condition  left  to  be  satisfied  in  order  that  M  be  a 
solution  of  the  Einstein  equation  for  a  T-wall  in  the  thin  shell 
approximation  is  equation  (1.10).  Using  equations  (2.18)  and 
(2.29)-(2.31)    yields 


Kab+"  Kab-=  vavb(d/dr){L+-  LJ  -  qabr"1{L+-  LJ  ,  (2.34) 


while  equations  (1.10)  and   (2.12)   give 

Kab+-  Kab-=  d/2)Ka{vavb(2r-l)  -  qab]  .  (2.35) 

Thus  Einstein's  equation  reduces  to  the  following  pair  of  equations: 
(d/dr){L+-  LJ  =  Ka[r-(l/2)]  ,  and  (2.36) 


r"1{L+-  LJ  =  (1/2)ko  .  (2.37) 


To   obtain   a  first  integral   of  these   equations,   multiply   equation   (2.36) 
by   Tl'2r,    multiply    equation    (2.37)    by    (l-2r)r1_2r,    and    add    the 
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resulting   two   equations   to   obtain 


r1"2r(d/dr){L+  -  L  }  +  (l-2r)r1_2r{L+-  LJ  =  0  ,  (2.38) 

(d/dr){r1"2r(L+-  L_)}  =  0  ,  and  (2.39) 

L+-  L.  =  c0r2r"1,  (2.40) 


where   cQ   is  a  constant.  Note  that  it  follows  from  equations  (2.37)  and 
(2.40)    that 


a  =  2c0K"1r2r-2.  (2.41) 


Since  the  walls  treated  here  have  positive  energy  density,  c  >  0.  Let 
F+  and  s+  be  the  values  of  F  and  s  on  I+,  with  F_  and  s_  defined 
similarly.  Using  equations  (2.23),  (2.25)  and  (2.40),  obtain 

s+(F+  +  (dRAft)2)1/2  -  s.(F.  +  (dRAR)2)1/2  =  c0R2rM.         (2.42) 

Now  using  the  general  interpretation  of  the  quantity  s,  it  is 
possible  to  interpret  the  quantities  s+  and  s_  as  follows:  since  na 
points  from  M"  to  M+,  the  unit  normal  pointing  into  M+  is  na  and  the 
unit    normal    pointing    into    M"  is   -na.  Hence  if  s ,  =   1    then  M+  is  an 
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interior;  if  s+=-l  then  M+  is  an  exterior.  Correspondingly  if  s_=  1  then 
M"  is  an  exterior;  if  s_=-l  then  M"  is  an  interior.  Since  cQ>  0,  it  follows 
from  (2.42)  that  the  case  s+=-l,  s_=   1   is  forbidden.  However,  all  other 

possible  combinations  are  allowed.  Thus  the  thin  wall  space-time  M 
can  consist  of  an  exterior  glued  to  an  interior  or  two  interiors  glued 
together,  but  cannot  consist  of  two  exteriors  glued  together. 

Next,   eliminate   s+  and  s_  from  the  equation  of  motion  for  R(X). 

Adding  s_(F_+  (dR/dX)2)1/2  to  both  sides  of  equation  (2.42)  and 
squaring    gives 

-S_(F.  +  (dRAR)2)1/2  = 
=  (1/2)(F_-  F+)c0-1R1-2r+(l/2)c0R2r-1.  (2.43) 

Squaring   equation   (2.43)  results   in 


(dR/d?02  =  {(F_-  F+)/2c0R2r"1}2  + 

+  (l/4)c02R4r'2  -  (1/2)(F++F_)  .  (2.44) 


A   solution   R(?i)   of  equation  (2.44)  determines   the  constants   s+  and  s. 
as  follows:  given  R(X,)  that  satisfies  (2.44),  solve  equation  (2.43)  for  s_. 
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Then  solve  (2.42)  for  s+.  The  function  R(X)  determines  the  surfaces  Z+ 
and  I".  The  quantity  s+  determines  which  region  of  M+  (the  exterior 
or  the  interior)  forms  with  Z+  the  manifold  M+.  The  corresponding 
result  holds  for  s_  and  M".  Note  that  s_  changes  sign  when  the 
right-hand  side  of  equation  (2.43)  changes  sign.  This  does  not 
correspond  to  any  pathological  behavior  of  the  surface  E"  or  of  the 
manifold  M",  but  is  simply  a  point  where  the  vectors  va  and  (3/3r)a 
in  M"  are  parallel.  The  corresponding  result  holds  for  s,.  Thus  given 
the  space-times  M+  and  M".  a  solution  of  equation  (2.44)  determines 
the  motion  of  the  T-wall. 

Using  equations  (2.2)  and  (2.44),  obtain 

(dRAR)2  =  {[(m+-  m_)  +  (l/6)(A+-A_)R3]/c0R2r}2  + 
+  (l/4)c02R4r"2-l  +  (m+  +  mJR"1    +  (1/6)(A+  +  A_)R2  .       (2.45) 

Define  the  constants  a,  b,  c,  f,  and  E  and  the  variable  X  by 


a  =    2(m+-  m_)/c02  ,  (2.46a) 

b=  (A  +  -A.)/3c02  ,  (2.46b) 

c=  4(m+  +  m.)/c02  ,  (2.46c) 

f  =  2(A+  +  A_)/3cQ2  ,  (2.46d) 
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E  =  -4/c02  ,  (2.46e) 

k  =  c0X/2  .  (2.46f) 


Then   equation   (2.45)   becomes 


(dR/cR)2  +  Y(R)  =  E  ,  (2.47) 


where  V_(R)  is  given  by 

V(R)  =  -(aR"2r+  bR3"2r)2  -  R4r"2-  cR"1  -  fR2  .  (2.48) 

Equation  (2.47)  is  just  the  equation  of  motion  for  a  particle  with 
energy  E  moving  in  one  dimension  in  a  potential  V_(R).  Here  R  plays 
the  role  of  the  spatial  coordinate,  L  plays  the  role  of  time,  and  the 
mass  of  the  "particle"  is  2.  Note  that  R  takes  on  only  positive  values 
and  that  E  takes  on  only  negative  values.  Thus  the  problem  of  the 
motion  of  the  T-wall  has  been  reduced  to  an  equivalent  problem  in 
one-dimensional  mechanics.  The  properties  of  the  one-dimensional 
particle  motion,  and  thus  of  the  T-wall  motion,  depend  on  the 
properties  of  the  function  V_(r)  and  on  the  value  of  the  constant  E.  Any 
local  minimum  of  V_(R)  represents  a  possible  point  of  stable 
equilibrium,  while  any  local  maximum  of  V_(R)  represents  a  possible 
point  of  unstable  equilibrium.   Any  point  where  V_(R)  =  E  represents  a 
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possible  turning   point  of  the   motion. 

The  properties  of  the  function  V_(R)  depend  on  the  values  of  the 
constants  a,  b,  c,  f,  and  I\  In  this  work,  certain  restrictions  have  been 
placed  on  the  values  of  these  constants.  Recall  that  0<  r<l,  and  that 
A+,  A.  ,  m+,  and  m_  are  all  non-negative.  It  then  follows  that 

2lal<  c  ,  and  (2.49a) 

2lbl<f.  (2.49b) 

Thus    the    constants    c    and    f    are    nonnegative.    It    then    follows    from 
(2.48)  that  V(R)<0  for  all  R. 

Various  special  cases  of  the  wall's  equation  of  motion  have  been 
treated  previously.  These  include  the  case  a=b=f=0  treated  in  5,7;  the 
case  r=l,  a=c=0  treated  in  8,  9,  10;  the  case  r  =  l,  a=c/2  treated  in 
11-16;  the  case  r=b=f=0,  a=c/2  treated  in  3;  and  the  case  r=l,  a=0 
treated  in  17.  In  addition,  a  numerical  treatment  of  the  general 
problem  has  been  done  in  18.  The  method  of  reference  18  is  to  choose 
various  values  of  the  parameters  (a,b,c,f,r)  and  numerically  produce 
plots  of  the  corresponding  potentials.  Since  the  potential  has  so  many 
parameters,  this  is  a  somewhat  unwieldy  way  to  obtain  information 
about  the  wall's  motion.  In  this  work,  analytic  methods  will  be  used 
to  determine,  as  a  function  of  the  parameters,  certain  features  of  the 
wall's  motion:  the  asymptotic  values  of  R,  whether  the  wall's  motion 
has    a    turning    point    and    whether    the    wall    oscillates    about    stable 
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equilibrium. 

The  possible  values  of  the  parameters  c  and  f  can  be  divided 
into  a  "generic"  case:  c*0,  f*0,  and  three  "exceptional"  cases:  (i)  c=0, 
f*0;  (ii)  c*0,  f=0,  and  (iii)  c=f=0. 

Case  1:  c*0.  f*0 

Under  these  conditions  it  follows  from  equation  (2.48)  that  V.  is 
unbounded  below  both  as  R-»0  and  as  R->».  This  means  V_(R)  must 
have  at  least  one  local  maximum.  In  fact,  it  is  possible  to  demonstrate 
that  there  is  a  number  Tc  such  that  if  T>TC  then  V_(R)  has  at  least  one 

local   maximum  and  no   other  stationary  points.   Assume  that     I>0.   Let 
Rs  be  a  stationary  point  of  V_(R),  and  let  a  prime  denote  differentiation 

with  respect  to  R.  Then     V_'(RS)  =  0.  It     follows  from  equation  (2.48) 
that 

V"(RS)  =  {R_1(d/dR)(RdVydR)}   = 
=  -16r2a2Rs"4r-2-(6-4r)2b2Rs4-4r-  2(3-4r)2abRs1"4r- 

-  cRs"3  -  4f  -  (4r-2)2Rs4r~4.  (2.50) 

Rearranging   the   terms   on   the  right-hand   side   of  (2.50)   results   in   the 
following    equation: 
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Y"(RS)  =  -i6r-2Rs"4r-2{ar2+  b(r-(3/4))2Rs3}2  + 
+  18r-2b2Rs4"4r{r  -    3(2-  /2)8_1 }  {r  -  3(2  +  /2)8_1 }  -  cRs"3  - 
-4f-(4r-  2)2Rs4r"4  .  (2.51) 

All  terms  on  the  right-hand  side  of  (2.51)  except  for  the  second  are 
manifestly  negative.  In  the  second  term,  since  T<1  it  follows  that  the 
factor  {r-  3(2+/ 2)8-1}  is  negative.  Now    define  Tc  by 

Tc=   3(2-/2)/8  .  (2.52) 

Then  if  r>rc  it  follows  that  the  second  term  on  the  right-hand  side  of 
equation  (2.51)  is  manifestly  negative.  Consider  the  case  T>rc  .  Then 
V"(RS)  is  negative.  Thus  Y.(RS)  is  a  local  maximum,  and  any  stationary 

point  of  V(R)  is  a  local  maximum.  Since  V_(R)  is  smooth  it  follows  that 
it  can  have  at  most  one  maximum,  and  since  it  has  at  least  one 
maximum,  it  has  exactly  one  maximum. 

The  one-dimensional  motion  of  a  particle  in  a  potential  with  one 
maximum  and  no  other  stationary  points  is  fairly  simple.  The  motion 
has  either  one  turning  point  or  no  turning  points  depending  on 
whether  the  energy  is  less  than  the  maximum.  Since  equation  (2.47) 
is    time-reversal    invariant,    for    any    possible    motion    of   the    wall,    the 
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time-reversed   motion   is   also   possible. 

Let  V_max  be  the  maximum  value  of  V_(R)  and  let  Rs  be  the  value 
of  R  which  maximizes  V_(R).  Then  the  possible  motions  of  the  wall 
depend  on  Ymax    and  E.  For  E>  V_max  the  possible  motions  of  the  wall 

are  as  follows:  (i)  the  wall  appears  with  zero  size  and  expands  forever 
to  unbounded  size;  (ii)  the  time-reverse  of  the  motion  in  (i).  For  E= 
—max  me  PosslDle  motions  of  the  wall  are  as  follows:  (i)  the  wall 
remains  in  unstable  equilibrium  with  R=R§;  (ii)  the  wall  appears  with 
zero  size  and  expands  with  R  asymptotically  approaching  R§;  (iii)  the 
wall  contracts  from  unbounded  size  with  R  asymptotically 
approaching     R§;    (iv)   the   time-reverse   of  the   motion   in   (ii);    (v)   the 

time-reverse  of  the  motion  in  (iii).  For  E<  V_max  the  possible  motions 
of  the  wall  are  as  follows:  (i)  the  wall  appears  with  zero  size,  expands 
to  a  maximum  R  less  than  Rs,  contracts  to  zero  size  and  disappears;  (ii) 
the  wall  contracts  from  unbounded  size  to  a  minimum  R  greater  than 
Rs  and  then  expands  to  unbounded  size. 

An    exhaustive   description    of   the   possible    motions    of   the   wall 
has  been  given,  in  the  generic  case,  for     T>  r*c.  The    T<  Tc   is  more 

difficult.  One  of  the  features  of  the  motion  in  the  r>  Tc  case  is  that 
stable  equilibrium  is  not  possible.  This  feature  is  not  always  shared 
by  the    T<  Tc  case.  That  is,  there  are  some  values  of  the  constants  a,  b, 
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c,  f,  E  and  T,  with    r<  Tc,  such  that  the  wall  is  in  stable  equilibrium.  In 

particular,  choose  r=0  and  let  b  be  greater  than  4.  Choose  a,c,  and  f  as 
follows:  a=-3b  ,  c=6b  ,  f=9b  .  Then  using  these  choices  in  equation 
(2.48)     results  in 

V'(R)  =  6b2(-R5+3bR2-3bR+R"2)  +  2R-3.  (2.53) 

It    follows  from  (2.53)  that      V'(l)=2,  so  that  in  particular 

V'(1)>0  .  (2.54) 

Now  define  a  by  aH(3b2)-1/8.  Note  that  cc<l.  Using  equation  (2.53), 
obtain 

V'(a)=6bV2(l+/3-35/8b1/4)  .  (2.55) 

Since  b>4  it  follows  that  V'(oc)<0.  This  means  V'  changes  from  negative 
to  positive  as  R  goes  from  a  to  1,  and  hence  there  must  be  a  local 
minimum  at  some  value  of  R  in  that  interval.  It  follows  that  for 
certain  energies  the  one-dimensional  motion  in  the  potential  V_(R)  is 
an  oscillation  about  a  point  of  stable  equilibrium.  The  wall,  therefore, 
can  undergo  spherical  ocillations  about  a  radius  of  stable  equilibrium. 
Since     V_(R)   is   smooth   in   the   parameters    (a,    b,   c,    f,  T),   any   set  of 
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parameters  sufficiently  near  the  ones  given  will  also  give  a  V_(R) 
which  has  a  local  minimum  and  thus  give  rise  to  walls  with  similar 
dynamics. 

Case  2:  c=0.  f*0 

Under  these  conditions,  it  follows  from  equation  (2.49a)  that 
a=0.  Thus  m+=m  _=0  and  the  space-times  M+  and  M"  are  both  regions 
of  de  Sitter  space-time.  From  equation  (2.48)  it  follows  that 

V(R)  =  -b2R6"4r  -  R4r"2  -  fR2.  (2.56) 

It  follows  from  (2.56)  that  V_(R)  is  unbounded  below  as  R-><».  It  also 
follows  from  (2.56)  that  the  behavior  of  V_(R)  as  R-»0  depends  on  the 
value  of  T.  For  T<l/2,  V(R)  is  unbounded  below  as  R->0.  For  r=l/2, 
limR^o  V(R)=-1.  For  r>l/2,  HmR_>0  V(R)=0.  For  T>l/2  it  follows  from 

equation  (2.56)  that  V(R)  has  no  stationary  points.  For  T<l/2  one  can 
show  using  (2.56)  that  V_(R)  has  exactly  one  maximum  and  no  other 
stationary  points.  In  this  latter  case,  let  Y.max  and  R§  be  respectively 
the  maximum  value  of  V(R)  and  the  value  of  R  at  which  V_(R)  attains 
this  maximum.  Then  the  possible  motions  of  the  wall  depend  on  T, 
—max'  anc*  E-  For  T<l/2,  E<  V_max  the  possible  motions  of  the  wall  are 
as  follows:   (i)  the  wall  appears  with  zero  size,  expands  to  a  maximum 
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R  less  than  Rg,  contracts  to  zero  size  and  disappears;  (ii)  the  wall 
contracts  from  unbounded  size  to  a  minimum  R  greater  than  R§  and 
then  expands  to  unbounded  size.  For  T<l/2,  E=  V_max  the  possible 
motions  of  the  wall  are:  (i)  the  wall  remains  in  unstable  equilibrium 
with    R=RS;   (ii)  the  wall   appears  with  zero  size  and  expands  with  R 

asymptotically  approaching  R§;  (iii)  the  wall  contracts  from 
unbounded  size  with  R  asymptotically  approaching  Rs;  (iv)  the 
time-reverse  of  the  motion  in  (ii);  (v)  the  time-reverse  of  the  motion 
in  (iii).  For  T<l/2,  E>V_max  the  possible  motions  of  the  wall  are  as 
follows:  (i)  the  wall  appears  with  zero  size  and  expands  to  unbounded 
size;  (ii)  the  time-reverse  of  the  motion  in  (i).  For  T  =  l/2,  E<-1  the 
motion  of  the  wall  is  as  follows:  the  wall  contracts  from  unbounded 
size  to  a  minimum  R  and  then  expands  to  unbounded  size.  For  T=l/2, 
E=-l  the  possible  motions  of  the  wall  are  as  follows:  (i)  the  wall 
contracts  from  unbounded  size  with  R  asymptotically  approaching 
zero;  (ii)  the  time-reverse  of  the  motion  in  (i).  For  r=l/2,  E>-1  the 
possible  motion  of  the  wall  is:  (i)  the  wall  contracts  from  unbounded 
size  to  zero  size  and  disappears;  (ii)  the  time-reverse  of  the  motion  in 
(i).  For  T>l/2  the  motion  of  the  wall  is  as  follows:  the  wall  contracts 
from  unbounded  size  to  a  minimum  R  and  then  expands  to 
unbounded    size. 
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Case  3.  c*0.  f=0 

From    equation    (2.49b)    it    follows    that    b=0.    Thus    A  +  =  A  =0; 

therefore  the  vacuum  energy  outside  the  wall  is  zero  and  the 
space-times  M+  and  M"  are  both  regions  of  Schwarzschild  space-time. 
It  then  follows  from  equation   (2.48)  that 

V(R)=  -a2R_4r-R4r"2-cR_1.  (2.57) 

It  follows  from  (2.57)  that  V(R)  is  unbounded  below  as  R->0.  For  r>l/2 
it  follows  from  equation  (2.57)  that  V_(R)  is  unbounded  from  below  as 
R->».  For  r=l/2,  lim^*,  V(R)=-1.  For  T<l/2,  lim^,,  V(R)=0.  For  T<l/2 
it  follows  from  equation  (2.57)  that  V_(R)  has  no  stationary  points. 
Once  can  show  using  equation  (2.57)  that  for  T>l/2,  V(R)  has  exactly 
one  maximum  and  no  other  stationary  points.  For  T>l/2  define  R§  and 

—max   as  before.  The  possible  motions  depend  on  r,  V_max,  and  E.  For 

T<l/2  the  motion  of  the  wall  is  as  follows:  the  wall  appears  with  zero 
size,  expands  to  a  maximum  R,  contracts  to  zero  size  and  disappears. 
For  T=l/2,  E<-1  the  motion  of  the  wall  is  as  follows:  the  wall  appears 
with  zero  size,  expands  to  a  maximum  R,  contracts  to  zero  size  and 
disappears.  For  T=l/2,  E  >  -1  the  possible  motions  of  the  wall  are  as 
follows:  (i)  the  wall  appears  with  zero  size  and  expands  to  unbounded 
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size;  (ii)  the  time-reverse  of  the  motion  in  (i).  For  r>l/2,  E<  V_._Y  the 


max 


possible  motions  of  the  wall  are  as  follows:  (i)  the  wall  appears  with 
zero  size,  expands  to  a  maximum  R  less  than  R§,  contracts  to  zero  size 
and  disappears:  (ii)  the  wall  contracts  from  unbounded  size  to  a 
minimum   R   greater  than   Rs   and  then  expands  to  unbounded   size.  For 

T>l/2,  E=V_max  the  possible  motions  of  the  wall  are  as  follows:  (i)  the 
wall  remains  in  unstable  equilibrium  with  R=RS;  (ii)  the  wall  appears 
with  zero  size  and  expands  with  R  asymptotically  approaching  R  ;  (iii) 
the  wall  contracts  from  unbounded  size  with  R  asymptotically 
approaching       Rg;   (iv)   the   time-reverse   of  the   motion   in   (ii);   (v)   the 

time-reverse  of  the  motion  in   (iii).  For  T>l/2,  E>  Y_max    the  possible 

motions  of  the  wall  are  as  follows:  (i)  the  wall  appears  with  zero  size 
and  expands  to  unbounded   size;   (ii)  the  time-reverse  of  the  motion  in 

(i). 

Case  4:  c=f=0 

From  equations  (2.49a)  and  (2.49b)  it  follows  that  a=b=0.  Thus 
A  +  =A  =m  +  =  m  =0  and  M+  and  M"  are  both  regions  of  Minkowski 
space-time.   It  then   follows   from  equation   (2.48)   that 


V(R)=-R4r"2.  (2.58) 
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For  T<l/2  it  follows  from  equation  (2.58)  that  lim^,,,  V(R)=0,  V(R)  is 
unbounded  from  below  as  R-+0  and  Y(R)  has  no  stationary  points.  For 
T=l/2,  it  follows  from  (2.58)  that  V(R)=-1.  For  r>l/2,  limR_,0  V(R)=0, 
V_(R)  is  unbounded  below  as  R-*co  and  \L(R)  has  no  stationary  points. 
The  dynamics  of  the  wall  depend  on  T  and  E.  For  T>l/2  the  motion  of 
the  wall  is  as  follows:  the  wall  contracts  from  unbounded  size  to  a 
minimum  radius  and  then  expands  to  unbounded  size.  For  T=l/2,  E=-l 
the  motion  of  the  wall  is  as  follows:  the  wall  remains  with  R  at  a  fixed 
value.  For  r=l/2,  E  >-l  the  possible  motions  of  the  wall  are  as  follows: 
(i)  the  wall  appears  with  zero  size  and  expands  to  unbounded  size;  (ii) 
the  time-reverse  of  the  motion  given  in  (i).  For  T<l/2  the  motion  of 
the  wall  is  as  follows:  the  wall  appears  with  zero  size,  expands  to  a 
maximum  radius,   contracts  to  zero  size  and  disappears. 

Global    Structure 

The  wall  space-time  consists  of  two  regions  of  Schwarzschild-de 
Sitter  space-time  attached  along  timelike  spherically  symmetric 
surfaces.  Thus  we  first  examine  the  global  structure  of 
Schwarzschild-de  Sitter  space-time.  The  expression  given  for  the 
Schwarzschild-de  Sitter  metric  in  equation  (2.1)  is  badly  behaved  at 
points  where  F=0;   however  the  curvature   scalar  Rahcd1*^0^    diverges 
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only  as  r-»  0.  Thus  the  coordinate  system  in  equation  (2.1)  covers  a 
region  of  the  manifold  where  F  is  nonzero  and  may  not  cover  the 
entire  space-time.  It  is  desirable  to  find  a  coordinate  system  which 
does  cover  the  entire  space-time.  To  do  this,  first  examine  the 
behavior  of  radial  null  geodesies. 

Consider   the    two-dimensional    metric 

ds2  =  -Fdt2  +  F_1dr2.  (2.59) 

Null  geodesies  in  this  metric  are  the  same  as  radial  null  geodesies  in 
the  metric  of  equation  (2.1).  In  the  two  dimensional  metric  of 
equation  (2.59)  let  ka  be  the  tangent  vector  to  a  null  geodesic  with 
affine  parameter     X.  Then  there  is  a  constant  co   such  that 


co=  ka(3/3t)a  =   -F(dtAR)  .  (2.60) 


Since  kak  =0  it  follows  that 

a 


0  =  -F(dtAR)2  +  F'l(dr/dX)2  = 
=  F"1{(dr/d?i)2-o)2}.  (2.61) 


Thus 


r  =  ±co^  +  cj,     and  (2.62a) 
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t  =  c2  ±  J-F_1dr  ,  (2.62b) 

where   Cj   and  c2  are  constants.  If  F  is  nonzero  for  all  values  of  r  then 

it  follows  from  equations  (2.62a-b)  that  for  any  incomplete  null 
geodesic,  r  approaches  zero  as  the  affine  parameter  approaches  the 
limiting  value.  If  F  is  zero  for  some  value  of  r,  let  rQ  be  the  smallest 
such  value.  The  it  follows  from  equations  (2.62a-b)  that  the 
space-time  has  incomplete  null   geodesies   for  which  r  approaches  r     as 

the   affine   parameter   approaches    the   limiting    value. 
Introduce  the  coordinates   u   and  v   by 

u  =  -t  +  /0r  F_1(r)dr,  (2.63a) 

w=  t  +  /0r  F^rjdr.  (2.63b) 

Then  the  metric  in  the  coordinates  (u,r)  is 

ds2  =  2dudr  -  Fdu2,  (2.64) 

and  in  the  coordinates  (v,r), 

ds2  =  2dvdr  -  Fdv2.  (2.65) 
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The  metric  in  equation  (2.64)  is  smooth  for  -  °°  <  u  <  °°  and  0<  r<  ». 
For  0<  r  <  rQ   the  metrics  in  equations  (2.64)  and   (2.59)  are  isometric. 

Thus  the  space-time  in  equation  (2.64)  is  an  extension  of  the 
space-time  in  equation  (2.59).  The  corresponding  results  hold  for  the 
metric  in  equation  (2.65).  Using  equations  (2.63a-b),  one  finds  that 
for  0  <  r  <  rQ  the  (u,r)  coordinate  patch  and  the  (v,r)  coordinate  patch 
overlap.  The  overlap  map  is 

v=-u  +  2/0rF-1(r)dr.  (2.66) 

This  larger  space-time  contains  both  the  (u,r)  coordinate  patch  and 
the  (v,r)  coordinate  patch.  In  general  it  takes  more  than  one  (u,r)  and 
one  (v,r)  coordinate  patch  to  cover  the  entire  space-time.  Introducing 
the  appropriate  (u,r)  and  (v,r)  coordinate  patches  and  appropriate 
overlap  maps,  it  is  possible  to  obtain  a  space-time  with  the  property 
that,  on  every  incomplete  null  geodesic,  r  approaches  0  as  the  affine 
parameter  approaches  its  limiting  value.  The  extended  space-time 
with   this  property   is  called   the  complete  extension. 

To  display  some  information  about  the  complete  extension  in  a 
simple  graphic  form,  introduce  coordinates  y  and  t,  such  that  the 
metric  has  the  form  ds2  =  Q2(-d£2  +  d\|/2)  for  some  function  Q  and 
such  that  all  null  geodesies  approach  finite  values  of  £  and  y  as  their 
affine   parameters   tend   to  ±°°    or  to  a  limiting  value.  Then  in  the  (£,\|/) 
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plane  draw  the  limit  points  of  the  null  geodesies  (including  those  at 
"infinite  affine  parameter").  These  points  can  be  regarded  as  a 
boundary  attached  to  the  space-time.  Also  draw  in  the  points  where 
F=0.  Points  representing  a  curvature  singularity  are  drawn  with  a 
sawtooth  line.  In  the  case  m=0  there  is  no  curvature  singularity  and 
the  points  at  r=0  are  drawn  with  a  dotted  line.  Each  point  of  the 
resulting  diagram  (called  a  Penrose  diagram)  can  be  regarded  as  a 
2-sphere  of  symmetry  in  Schwarzschild-de  Sitter  space-time.  Since 
radial  null  geodesies  are  lines  which  make  an  angle  of  45°  with  the 
vertical,  it  follows  that  the  path  of  the  wall,  which  is  a  timelike  curve 
on  this  diagram,  always  makes  an  angle  of  less  than  45°  with  the 
vertical.  In  drawing  the  Penrose  diagram,  choose  a  time  orientation 
for  the  space-time.  Choosing  the  opposite  time  orientation  simply 
results  in   the  same  diagram   "upside  down." 

The  nature  of  the  complete  extension  depends  on  the  values  of 
m  and  A.  Recall  that  only  the  cases  where  A  and  m  are  nonnegative 
are  considered  here.  There  are  two  "generic"  cases:  (i)  9m  A  >  1,  (ii) 
0<9mz,A<l;  and  four  "exceptional"  cases:  (i)  9irrA  =  l,  (ii)  m=A=0,  (iii) 
m=0,  A*0,  (iv)  m*0,  A=0.  The  second  generic  case  is  covered  in 
reference  25  and  last  three  exceptional  cases  are  covered  in  reference 
26,  so  the  results  for  these  cases  will  be  only  briefly  presented.  Some 
information  about  the  first  generic  and  first  exceptional  cases  is 
presented    in    reference    16. 

In  the  generic  case  9m2A>l,   it  follows  from  equation   (2.2)  that 
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F<0    for    all    r.    Thus    the    (r,t)    coordinate    patch    covers    the    complete 
extension.   Define  the  coordinates  %  and  y  and  the  constant    £„     by 

Xf  =  t  (2.67) 

%  s  -VF-^dr  (2.68) 

U  s  limr_.t  (2.69) 

Then  it  follows  from  equation  (2.59)  that 

ds2  =  -F(-d^2  +  dy2).  (2.70) 

The  points  at  b,=0  correspond  to  r=0;  those  at  ^=^o>    correspond  to  r=°° . 

The  Penrose  diagram  of  this  space-time  is  shown  in  figure   1,  page  42. 
The  lines  are  labelled  with  the  corresponding  values  of  r. 

In  the  exceptional  case  9m2A=l,  it     follows  from  equation  (2.2) 
that  F  and   (dF/dr)   vanish   at  r=rQ   where  rQ=3m.   For  all  r*rQ,    F(r)<0. 

The  complete  extension  is  covered  with  coordinate  patches  (u_,  r)  and 
(vn,  r),  one  for  each  integer  n.  The  metric  on  each  (un,  r)  coordinate 
patch  has  the  form  given  in  equation  (2.64)  and  in  each  (vn,  r) 
coordinate  patch  the  metric  has  the  form  given  in  equation  (2.65). 
Each    (un,    r)    coordinate   patch    overlaps    with    the    (vn,    r)    coordinate 
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patch   for  r>rQ  and  with  the  (vn  +  1,  r)  coordinate  patch  for  r<  r  .   The 
overlap  maps  are,  for  r>rQ 


•   rrl 


vn  =  -un-2/r     F-1(r)dr,  (2.71) 


and  for  r<  r 


vn+1  =-un  +2/0rF-1(r)dr.  (2.72) 

Using  the  geodesic  equation  one  can  show  that  all  null  geodesies  in 
this  space-time  either  have  r=rQ  and  are  complete  or  have  r  take  on 
all  values  between  0  and  °°   and  are  incomplete  only  as  r->0. 

Next,    introduce   the    coordinates    £,  and  y   as  follows:  for  0<  t,<  k 
define       un    and      vn   in     terms     of     £     and    \\t       as      follows:      for 

(27tn-£)<\|/<27t(n+lK 

un=  cot{(^+V)/2}.  (2.73) 

For  (2ji(n-l)  +  ^}<\|/<(27in+^) 

vn=  cot{(£-v)/2}.  (2.74) 
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Define    r  in  terms  of  %  and  \|/  as  follows:  for  \j/  +  l=lTin  or  \j/-£=27in,  r=r 

For  other  values  of  \|/  and  £  (with  0<£<7t)  define  r  by  equations  (2.73), 
(2.74)   and  either  (2.71)   or  (2.72).   Then   it  follows   that  un  and  vn   are 

C°°  functions  of  £  and  \\f  and  that  r  is  a  C1  function  of  ^  and  \\f .  The 
metric  in  the  E,    and  y   coordinates  is 

ds2=  -F(-d^2 +dy2){4sin2[(^+V)/2]sin2[(^-y)/2]  J'1,       (2.75) 

where  r  is  regarded  as  a  function  of  I,  and  \\r  and  the  metric 
coefficients  are  defined  by  continuity  at  those  points  where  r=rQ.  The 
Penrose  diagram  for  this  space-time  is  shown  in  figure  2,  page  43. 
The  geodesic  limit  points  with  r=0  occur  at  £=0,  \j/*27tn;  those  with  r=°° 
occur  at  t=%,  \j/*(2n  +  l)7i;  those  with  r=rQ  occur  at  ^=0,  \jr=2rcn  and  ^=k, 
y=(2n+l)7i. 

The  Penrose  diagram  for  the  space-time  with  0<9m2A<l  is 
shown  in  figure  3,  page  44.  5  In  this  case,  F(r)  vanishes  at  two  values 
of  r,  designated     r+  and  r_. 

Finally  the  results  for  the  last  three  exceptional  cases  will  be 
presented.  °  The  case  m=A=0  is  Minkowski  space-time.  Here  F=l  and 
there  is  no  curvature  singularity.  The  Penrose  diagram  for  Minkowski 
space-time   is   shown   in   figure   4,   page   45.   The   case   m=0,   A*0  is  de 
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Sitter  space-time,  with  F=0  at  r=(3/A)ly'2  and  no  curvature 
singularity  (figure  5,  page  46).  The  case  m*0,  A=0  is  Schwarzschild 
space-time,  where  F=0  at  r=2m  (figure  6,  page  47). 

There  are  many  possibilities  for  the  global  structure  of  the  wall 
space-time.  Each  of  the  manifolds  M+  and  M"  is  a  region  of 
Schwarzschild-de  Sitter  space-time  with  a  boundary  which  is  a 
timelike  spherically  symmetric  surface.  Thus  the  Penrose  diagram  of 
M  or  M"  consists  of  a  piece  of  one  of  the  Penrose  diagrams  of  figures 
1-6  bounded  by  a  timelike  curve.  The  Penrose  diagram  of  the  wall 
space-time  consists  of  the  Penrose  diagrams  for  M+  and  M"  with  the 
boundaries  identified  at  corresponding  values  of  r.  Figure  7,  page  48, 
is   the  Penrose  diagram  for   one  possible   wall   space-time.   In   this   case 

9m+  A  +  >1  and  0<9m_  A_<1.  The  curved  dotted  lines  represent  the 
paths  of  the  wall  in  M+    and  M"    and  are  identified. 
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r=oo 


r=0 


Figure   1.   The  Penrose  diagram  for  Schwarzschild-de   Sitter 
space-time  in  the  case  9irrA>l. 
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r=oo 


r=0 


r=0 


Figure  2.   The  Penrose  diagram  for  Schwarzschild-de   Sitter 
space-time  in  the  case  9m2A=l. 
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r=0 
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r=oo 


r=oo 


r=oo 
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r=o 


Figure  3.   The  Penrose  diagram  for  Schwarzschild-de   Sitter 
space-time  in  the  case  0<9m2A<l. 


4  5 


t=oo 


r=oo 


r=0 


Figure  4.  The  Penrose  diagram  for  Minkowski   space-time 
(m=  A  =  0). 
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r=0 
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r=oo 


Figure  5.  The  Penrose  diagram  for  de  Sitter  space-time 
(m=0,  A*  0).  Here,  rQ=(3/A)1/2. 
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r=0 


Figure  6.   The  Penrose  diagram  for  Schwarzschild   space-time 
(m*0,  A=0).  r  =  2m  and  t  =  *   on  the  inner  diagonals. 
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r=0 


r=0 


Figure  7.  The  Penrose  diagram  for  a  possible  spherical 
wall  space-time.  9m+2A>l,  0<9m_  A<1.  The  dotted  lines 

are    identified. 
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Planar  Walls 

To  find  the  motion  of  a  planar  wall  one  must  first  find  two 
solutions  of  the  vacuum  Einstein  equation  with  cosmological  constant 
with  planar  symmetry.  A  space-time  has  planar  symmetry  if  its 
symmetry  group  contains  the  symmetry  group  of  the  Euclidean 
2-plane.  The  general  planar  symmetric  solution  of  the  vacuum  Einstein 
equation   with   cosmological   constant  A  has  the  form 

ds2  =  -Hdp2  +  2drdp  +  r^dx2  +  dy2)  ,  (2.76) 

where  the  function  H  is  defined  by 

H  =    y/r  -   Ar2/3  (2.77) 

and  y  is  a  constant.  The  surfaces  of  constant  r  and  (3  are  surfaces  of 
planar  symmetry.  In  the  case  y*0  this  space-time  has  a  curvature 
singularity  at  r=0,  while  y=0  corresponds  to  de  Sitter  space-time,  and 
Y=A=0  is  Minkowski  space-time. 

There  is  another  type  of  planar  symmetry  in  Minkowski 
space-time.  Writing  the  Minkowski  metric  in  the  usual  Cartesian 
coordinates, 


ds2  =  -dt2  +  dz2  +  dx2  +  dy2,  (2.78) 
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it  is  manifest  that  the  surfaces  of  constant  t  and  z  are  surfaces  of 
planar  symmetry.  However,  the  volume  element  on  these  surfaces  is 
independent  of  t  and  z.  This  form  of  planar  symmetry  is  therefore  not 
suitable  for  describing  an  expanding  or  contracting  wall  and  will  not 
be  considered      here. 

In  the  case  where  y*  0  or  A  *  0  one  can  rewrite  the  metric  in 
equation  (2.76)  by  introducing  the  coordinate  p  defined  by 
p=~P+/H"  dr.  Then  the  metric  in  equation  (2.76)  is  given  by 

ds2  =  -Hdp2  +  H_1dr2  +  r^dx2  +  dy2)  .  (2.79) 

This  is  similar  to  the  expression  in  equation  (2.1)  for  the  spherically 
symmetric    case. 

Choose  M_  and  M_  to  be  space-times  with  metrics  of  the  form 
given  in  equation  (4.1).  Let  A+  and  y+  be  the  values  of  the  constants  A 

and    y  in  M.+  ,  respectively,  and  in  the  same  way  define  A_  and  y _.   Here, 

only  the  cases  where  both  A+  and  A_  are  greater  than  or  equal  to  zero 
will  be  considered. 

The  procedure  for  treating  planar  walls  and  deriving  their 
equation  of  motion  is  completely  analogous  to  the  method  used  in  the 
last  section  for  spherical  walls.  Therefore  we  omit  the  derivation  and 
simply   present   the   results.    The    wall    hypersurface    is    specified    by    a 
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relation  of  the  form 

T  =  R(k),  (2.80) 

where  X  is  proper  time  on  a  world-line  tangent  to  the  wall  and 
orthogonal  to  the  2-planes  of  symmetry.  The  energy  density  of  the 
wall  is  given  by 


a  =  2cn  K_1R2r-2  ,  (2.81) 


where  cQ  is  a  constant.  The  equation  of  motion  of  the  wall  is 

s+{H+  +  (dRAR)2  }1/2  -  s.{H.  +  (dR/<R)2}1/2=  c0R2rM.    (2.82) 

The  quantities  H+  and  H  denote  the  values  of  H  on  E+  and  Z" 
respectively.  The  constants  s+  and  s_  are  defined  as  in  the  spherical 
case.  From  equation  (2.82)     obtain 

(dRAR)2  =  (H_-  H+)2(2c0R2rMr2  +  c02R4r"2/4     - 

-(H+  +  H_)/2  .  (2.83) 

Then,  using  equation  (2.77),     obtain 


52 


(dR/dA.)2-  {(y_- Y+)+(l/3)(A+-  A.)R3}2(2c0R2r)-2  +  (l/4)c02R4r-2   - 
-(1/2)(Y++Y.)R"1    +  (1/6)  (A+ +  A.)R2  .  (2.84) 


Define  the  constants  a,  b,  c,  and  f  and  the  variable  ^  by 

a  =   (y.-Y+)/c02,  (2.85a) 

b=  (A+-  A.)/3c02,  (2.85b) 

c^  -2(T++Y.)/c02,  (2.85c) 

fs  2(A++ A.)/3c02,  (2.85d) 

k=  c0X/2  .  (2.85e) 

Then   equation   (2.84)   becomes 

(dRAR)2  +  V(R)  =  0  ,  (2.86) 
where  V(R)  is  given  by 

V(R)  s  -{aR"2r+  bR3_2r)2  -  R4r"2  -  cR1  -  fR2.  (2.87) 

Equation    (2.86)    is    the    equation    of    motion    for    a    particle  with    zero 
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energy  moving  in  one  dimension  in  a  potential  V_(R).  The  variable  R 
takes  on  values  in  the  open  interval  (0,  •).  The  special  case  a=b=f=0 
has  been  treated  in  references  4,5,6,  and  7. 

Since  A+  and  A_  are  nonnegative,  it  follows  that 

2lbl<f.  (2.88) 

Thus  the  constant  f  is  nonnegative  also.  Note,  however,  that  there  is  no 
restriction  on  the  sign  of  the  constant  c.  Consider  the  case  c  >0.  It  then 
follows  from  equation  (2.87)  that  V1(R)<0  for  all  R.  The  motion  of  a 
particle  of  zero  energy  has  no  turning  points  in  a  potential  that  is 
everywhere  negative.  If  a,  c,  and  r  satisfy  the  special  conditions  a=c=0, 
r>l/2,  then  lim^^Q  V_(R)=0  and  the  possible  motions  of  the  wall  are  as 

follows:  (i)  the  wall  contracts  from  unbounded  values  of  R  with  R 
asymptotically  approaching  zero;  (ii)  the  time  reverse  of  the  motion  in 
(i).  If  these  special  conditions  on  a,  c,  and  r  are  not  satisfied  then  the 
possible  motions  of  the  wall  are  as  follows:  (i)  the  wall  appears  with 
R=0  and  expands  forever  to  unbounded  values  of  R;  (ii)  the  time 
reverse  of  the  motion  in  (i). 

Next  consider  the  case  c<0.  Introduce  the  variable  Z  by 


Z=  R3/2.  (2.89) 
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Note  that  Z  is  a  monotonic  function  of  R;  thus  dZ/dL  has  the  same  sign 
as    dR/dA-   Rewriting  equation  (2.86)  in   terms   of  the  variable  Z  results 
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(4/9)(dZ/dl)2  +  W(Z)  =  c  ,  (2.90) 

where  the  function  W(Z)  is  defined  by 

W(Z)s  -z(2"8r)/3{a+bZ2}2  -z(8r"2)/3-  fZ2.  (2.91) 

Equation  (2.90)  is  the  equation  of  motion  of  a  particle  of  mass  8/9  and 
energy  c  moving  in  a  potential  W(Z). 

Next,  examine  the  properties  of  the  function  W(Z).  Consider  a 
"generic"  case:  r*l/4,  a*0  and  three  "exceptional"  cases:  (i)  r=l/4,  a*0; 
(ii)  r*l/4,  a=0  and  (iii)  r=l/4,   a=0. 

Case  1:  r*l/4.  a*0 

From  equation  (2.91)  it  follows  that  W(Z)  is  unbounded  from 
below  both  as  Z->0  and  as  Z->°°.  Thus  W(Z)  must  have  at  least  one  local 

maximum.  There  is  a  number  Tc  such  that  if  r>Tc  then  W(Z)  has 
exactly  one  local  maximum  and  no  other  stationary  points.  To  see  this, 
let    Zg     be    a    stationary    point    of    W(Z)    and    let    a    prime    denote 
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differentiation   with  respect  to  Z.   Then   W'(Zs)=0.   It  then  follows  from 


equation    (2.91)    that 


W"(ZS)={Z-1[ZW']'}(ZS)  = 
=  -(4/9)  Zs(-4-8r)/3(a2(4r_1)2+2ab(4r_4)2Zs2+  b2(4r-7)2  Z§4}  - 

-(4/9)(4r-l)2  Z  (8r-8)/3  .  4f  (2  92) 


Rearranging   terms   on   the  right-hand   side,      arrive  at 

W"(  Zs)  =  -(4/9)(4r-l)2Zs(-4-8r)/3{a+b(4r_4)2(4r.1)-2  ^2^2  + 

+  128b2(4r-l)"2  Zs(8-8r)/3{r.(8+3/2)8-1}{r-(8-3/2)8-1}- 

-(4/9)(4r-l)2  zs(8r'8)/3  -  4f  .  (2.93) 

All  terms  on  the  right-hand  side  of  (2.93)  are  manifestly  negative 
except  for  the  second  term.  In  the  second  term,  since  r<l  it  follows 
that  the  factor  r-(8+3/2)8"1   is  negative.  Now  define  the  number  Tc  by 

rc=  (8-3/2)/8  .  (2.94) 

Then  if  T>rc  it  follows  that  the  second  term  on  the  right-hand  side  of 
(2.93)  is  manifestly  negative.   Consider  the  case  r>rc.  Then  W"(  Z§)  is 
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negative.   Thus  W(  Zs)  is  a  local  maximum.  This  means  any  stationary 

point  of  W(Z)  is  a  local  maximum.  Since  W(Z)  is  smooth  it  follows  that 
it  can  have  at  most  one  maximum,  and  since  it  has  at  least  one 
maximum,   it  has  exactly  one  maximum. 

Let    Wmax    be  the   maximum  value  of  W(Z)   and   let  Zs   be  the 

value  of  Z  at  which  W(Z)  attains  this  maximum.  Define  R§  by  R§=  Z§2^. 
For  c>Wmax  the  possible  motions  of  the  wall  are  as  follows:  (i)  the 
wall  appears  with  R=0  and  expands  to  unbounded  values  of  R;  (ii)  the 
time  reverse  of  the  motion  in  (i).  For  c=  Wmax  the  possible  motions  of 
the  wall  are  as  follows:  (i)  the  wall  is  in  unstable  equilibrium  with 
R=RS;  (ii)  the  wall  appears  with  R=0  and  expands  with  R  asymptotically 

approaching  Rs;  (iii)  the  wall  contracts  from  unbounded  values  of  R 
with  R  asymptotically  approaching  Rs;  (iv)  the  time  reverse  of  the 
motion  in  (ii);  (v)  the  time  reverse  of  the  motion  in  (iii).  For  c<  Wmax 
the  possible  motions  of  the  wall  are  as  follows:  (i)  the  wall  appears 
with  R=0,  expands  to  a  maximum  R  less  than  R§,  contracts  to  R=0  and 
disappears;  (ii)  the  wall  contracts  from  unbounded  values  of  R  to  a 
minimum  R  greater  than  R§  and  then  expands  to  unbounded  values  of 
R 

The    generic    case    when    r<T      is   more  difficult,   but,   as   in   the 
spherical   case,   it   is   possible   to   demonstrate   that   for   some   choices   of 
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the  parameters  a,  b,  c,  f,  and  T  the  wall  can  undergo  oscillations  about 
a  point  of  stable  equilibrium.  Let  T  be  any  number  satisfying  0<T<l/4. 
Choose  a=-2b  and  f=b  .  Note  that  for  this  choice  of  f  to  be  compatible 
with  equation  (2.88)  one  must  have  b2>4.  Then,  using  equation  (2.81), 
find 

W'(Z)=-(2/3)b2Z-(1+8r)/3(Z2_2){(7_4r)Z2+  (8r-2)}- 

-2b2Z  +  (2/3)(l-4r)z(8r"5)/3  .  (2.95) 

From  (2.95)  it  follows  that 


W'(l)  =  (2/3)b2(2+4f)  +  (2/3)(l-4f),  (2.96) 


so  W'(1)>0.  Define  Z{  by  Zf{(2-Sr)/(7-4r)}l/2.  Then  Z,<1  and 

W,(Z1)=-2Z1{b2-(l/3)(l-4r)z/8r-8)/3}.  (2.97) 

Choose  b  so  that  b2>(l/3)(l-4r)  Z^Sr'S^3.  Then  from  (2.97)  it  follows 
that  W'(  Zj)  <  0.  Therefore  it  follows  that  W(Z)  has  a  local  minimum  at 
some  Z  between  Zj  and  1.  Hence  for  some  values  of  the  "energy"  c  and 
initial  value  of  R  the  wall  undergoes  oscillations  about  a  point  of  stable 


equilibrium. 

Case  2:  r=l/4.  a*0 

From  equation  (2.91)  it  follows  that 

W(Z)  =  -(a+bZ2)2  -  1  -  fZ2.  (2.98) 

From  equation  (2.98)  the  following  properties  hold  for  the  function 
W(Z):  (i)  if  (2ab+f)<0  then  W(Z)  has  exactly  one  local  maximum  and  no 
other  stationary  points;  (ii)  if  (2ab+f)>0  then  W(Z)  has  neither  local 
maxima  nor  local   minima;   (iii)   limz_>0W(Z)  =  -(l+a2);  (iv)  if  f=0  then 

(from  equation  (2.88))  b=0  and  W(Z)=-(l+a2).  Recall  that  Z  takes  on 
values  only  in  the  open  interval  (0,  •),  so  one  must  take  limits  as  Z-»0 
and  Z=0  cannot  considered  as  a  possible  location  for  a  minimum  or 
maximum  of  W(Z).  For  (2ab+f)<0  let  Wmax   be  the  maximum  value  of 

W(Z),  let  Z§  be  the  value  of  Z  at  which  this  maximum  is  achieved,  and 

let  Rs  be  given  by  RS=ZS2/3.  The  dynamics  of  the  wall  depend  on  2ab+f 

and  c.  For  (2ab+f)<0,  c>Wmax   the  possible  motions  of  the  wall  are  as 

follows:  (i)  the  wall  appears  with  R=0  and  expands  to  unbounded 
values  of  R;  (ii)  the  time-reverse  of  the  motion  in  (i).  For  (2ab+f)<0, 
c=Wmax  the  possible  motions  of  the  wall  are  as  follows:  (i)  the  wall  is 
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in  unstable  equilibrium  with  R=RS;  (ii)  the  wall  appears  with  R=0  and 
expands  with  R  asymptotically  approaching  R§;  (iii)  the  wall  contracts 
from  unbounded  values  of  R  with  R  asymptotically  approaching  R  ; 
(iv)  the  time-reverse  of  the  motion  in  (ii);  (v)  the  time-reverse  of  the 
motion  in  (iii).  For  (2ab+f)<0,  -(l+a2)<  c  <Wmax  the  possible  motions  of 
the  wall  are  as  follows:  (i)  the  wall  appears  with  R=0,  expands  to  a 
maximum  value  of  R  less  than  Rs,  contracts  to  R=0  and  disappears;  (ii) 
the  wall  contracts  from  unbounded  values  of  R  to  a  minimum  value  of 
R  greater  than  R§  and  then  expands  to  unbounded  values  of  R.  For 
(2ab+f)  <0,  c  <-(l+a  )  the  wall  contracts  from  unbounded  values  of  R 
to    a    minimum    value    of    R    greater    than    Rs     and    then    expands    to 

unbounded  values  of  R.  For  (2ab+f)  >0,  c  >-(l+a2)  the  possible  motions 
of  the  wall  are:  (i)  the  wall  appears  with  R=0  and  expands  forever  to 
unbounded  values  of  R;  (ii)  the  time-reverse  of  the  motion  in  (i).  For 
(2ab+f)>0,  c<-(l+a  )  the  wall  contracts  from  unbounded  values  of  R  to  a 
minimum  value   of  R   greater  than   R§    and   then   expands   to   unbounded 

values  of  R.  For  (2ab+f)>0,  f*0,  c=-(l+a2)  the  possible  motions  of  the 
wall  are  as  follows:  (i)  the  wall  contracts  from  unbounded  values  of  R 
with  R  asymptotically  approaching  zero;  (ii)  the  time-reverse  of  the 
motion  in  (i).  For  f=0,  c=-(l+a2)  the  wall  remains  with  R  at  any 
constant    value. 
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Case  3:  r*l/4.  a=0 

From  equation  (2.91)  it  follows  that 

W(Z)  =  -b2z(14-8r)/3  -  Z(8r-2)/3  -  fZ2.  (2.99) 

From  (2.99),  the  following  properties  hold  for  the  function  W(Z):  (i)  if 
T>l/4  then  limz_>QW(Z)=0  and  W(Z)  has  no  stationary  points;  (ii)  if 
T<l/4  and  f=0  then  b=0,  lim^.  W(Z)=0  and  W(Z)  has  no  stationary 
points;  (iii)  if  T<l/4  and  f*0  then  W(Z)  has  exactly  one  maximum  and 
no  other  stationary  points.  For  T<l/4,  f*0  define  Wmax,  Z  ,  and  R    as 

before.  The  dynamics  of  the  wall  depend  on  the  values  of  I\  f,  and  c. 
For  T>l/4  the  wall  contracts  from  unbounded  values  of  R  to  a 
minimum  value  of  R  and  then  expands  to  unbounded  values  of  R.  For 
T<l/4,  f=0  the  wall  appears  with  R=0,  expands  to  a  maximum  value  of 
R,    contracts    to    R=0    and    disappears.    For   r<l/4,    f*0,    c>Wmax    the 

possible  motions  of  the  wall  are:  (i)  the  wall  appears  with  R=0  and 
expands  forever  to  unbounded  values  of  R;  (ii)  the  time-reverse  of  the 
motion  in  (i).  For  T<l/4,  f*0,  c=Wmax  the  possible  motions  of  the  wall 
are:(i)  the  wall  remains  in  unstable  equilibrium  with  R=RS;  (ii)  the  wall 
appears   with   R=0   and   expands   with   R   asymptotically   approaching   R§; 
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(iii)  the  wall  contracts  from  unbounded  values  of  R  with  R 
asymptotically    approaching    R§;  (iv)  the  time-reverse  of  the  motion  in 

(ii);  (v)  the  time-reverse  of  the  motion  in  (iii).  For  T<l/4,  f*0,  c<W 
the  possible  motions  of  the  wall  are:  (i)  the  wall  appears  with  R=0, 
expands  to  a  maximum  value  of  R  less  than  Rs,  contracts  to  R=0  and 
disappears;  (ii)  the  wall  contracts  from  unbounded  values  of  R  to  a 
minimum  value  of  R  greater  than  R§  and  then  expands  to  unbounded 
values  of  of  R. 

Case  4:  r=l/4.  a=0 

From  equation  (2.91)  it  follows  that 

W(Z)=-b2Z4  -  1  -fZ2  .  (2.100) 

From  (2.100)  the  following  properties  hold  for  the  function  W(Z):  (i) 
limz_>0W(Z)=-l;  (ii)  if  f=0  then  W(Z)=-1;  (iii)  if  f*0  then  W(Z)  has  no 

stationary  points.  The  dynamics  of  the  wall  depend  on  the  values  of  f 
and  c.  For  c>-l  the  possible  motions  of  the  wall  are:  (i)  the  wall 
appears  with  R=0  and  expands  forever  to  unbounded  values  of  R;  (ii) 
the  time-reverse  of  the  motion  in  (i).  For  c=-l,  f *  0  the  possible 
motions  of  the  wall  are:  (i)  the  wall  contracts  from  unbounded  values 
of  R   with   R   asymptotically   approaching   zero;    (ii)   the   time-reverse   of 
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the  motion  in  (i).  For  c=-l,  f=0  the  wall  remain  at  any  constant  value  of 
R.  For  c<-l,  f  must  be  nonzero  and  the  wall  contracts  from  unbounded 
values  of  R  to  a  minimum  value  of  R  and  then  expands  to  unbounded 
values  of  R. 

Global    Structure 

The  global  structure  of  the  wall  space-time  consists  of  two 
regions  of  space-times  of  the  form  given  in  equation  (2.76)  attached 
along  timelike  plane  symmetric  surfaces.  Begin  by  considering  the 
properties  of  plane  orthogonal  null  geodesies,  that  is,  null  geodesies 
which  are  orthogonal  to  the  symmetry  2-planes.  Consider  the  two 
dimensional    metric 


ds2  =  2drdp  -  Hdp2.  (2.101) 


Null  geodesies  in  this  metric  are  the  same  as  plane-orthogonal  null 
geodesies  in  the  metric  of  equation  (2.76).  If  H  is  nonzero  for  all  values 
of  r  then  one  can  show  that  for  any  incomplete  null  geodesic  in  the 
metric  of  equation  (2.101)  r  approaches  zero  as  the  affine  parameter 
approaches  its  limiting  value.  If  H  vanishes  at  some  value  of  r, 
introduce  new  coordinates  to  extend  the  space-time.  In  either  case, 
obtain  a  space-time,  the  complete  extension,  in  which  any  incomplete 
null    geodesic    has    the   property    that   r   approaches    zero    as    the    affine 
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parameter   approaches   its    limiting    value. 

To  illustrate  the  properties  of  the  complete  extension,     introduce 

coordinates        £,     and    \|/      such     that     the     metric     has     the     form 

2       2         2         2 
ds   =Q   (-d£   +d\|/   )  for  some  function  Q   and  all  null  geodesies  approach 

finite  values  of  £  and  \y   as  their  affine  parameters  tend  to  ±°°     or  to  a 

limiting    value.    Then    produce    the    Penrose    diagram    of   the    space-time 

by  drawing,  in  the  £,\j/  plane,  the  limit  points  of  the  null  geodesies  and 

the  places  where  H  vanishes.  Each  point  in  the  Penrose  diagram  can  be 

regarded    as    a    2-plane    in    the    four-dimensional    space-time.    Curvature 

singularities   are   again  represented   by   sawtooth   lines;   if  y=0   then   there 

is   no   curvature    singularity   and   the   places   where   r=0   are   represented 

by   broken   lines. 

The  nature  of  the  complete  extension  depends  on  the  values  of  y 
and  A.  There  are  two  "generic"  cases:  (i)  y<0>A>0;  (ii)  Y>0>  A>0;  and  four 
"exceptional"  cases:  (i)  y=A=0;  (ii)  y=0,  A>0;  (iii)  y>0,  A=0;  (iv)  y<0,  A=0. 

In  the  generic  case  y<0,  A>0,  it  follows  from  equation  (2.77)  that 
H<0  for  all  r.  Thus  the  (P,r)  coordinate  patch  covers  the  entire 
space-time.   Define  the  coordinates  t,  and  \j/  and  the  constant  |„  by 

%  s   -Jo'H^Ctfdi,  (2.102) 

VE  P  +  ^  ,  (2.103) 

Us  limr^«^.  (2.104) 
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It  then  follows  from  equation  (2.101)  that 

ds2  =  -H(-d£2  +d\j/2)  .  (2.105) 

The  points  at  ^=0  correspond  to  r=0;  those  at  ^=^„  correspond  to  r=  °° 
The  Penrose  diagram  of  this  space-time  is  shown  in  figure  8,  page  70. 
In  the  generic  case  y>0,  A>0,  define  the  constant  Cj   by 

Cj  =  {3y/A}1/3.  (2.106) 

Then  it  follows  from  equation  (2.77)  that 


H=  (A/3r)(Cl3  -r3).  (2.107) 


Thus  H(cj)=0.  For  r<  cj,  H>0;  for  r>cl5  H<0.  The  coordinates  in  equation 

(2.101)  only  cover  part  of  the  space-time,  so     new  coordinates  must  be 
introduced  to  extend  the  space-time.  Define  the  function   S(r)  by 

S(r)  =  (4/3ACl2r)(r2+  Cjr  +  Cj2)3/2  X 

X  exp{/3{(jt/6)-tan-1[(2r+c1)/c1/3]}}  .  (2.108) 


65 
Then  for  r<  C|    it  follows  from  equations  (2.107)  and  (2.108)  that 

(d/dr){ln(S/H)}  =Ac1H"1.  (2.109) 

For  r<  Cj   define  the  coordinates  A  and  B  by 

A=  expf-Acjp/2},  (2.110) 

AB=  (2/ACl)2H/S  .  (2.111) 

Then,   using   equations   (2.101)   and   (2.109-2.111),      obtain 

ds2  =  SdAdB  .  (2.112) 

Here  r  is  a  function  of  A  and  B   given  by  equation  (2.111).  The  metric 
in  (2.112)  is  smooth  for  all  values  of  A  and  B.  The  region  A>0,  B>0  is 
isometric  to  the  (p,  r)  coordinate  patch.  This  space-time  is  the  complete 
extension  of  the  space-time  of  metric   (2.101). 
Define  the  coordinates  £,  and  \j/  by 

A  =  tan{(\|/-£)/2}  ,  (2.113) 

B  =  tan{(y+^)/2}   .  (2.114) 
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Then,   using  equations   (2.111)  and   (2.112),   obtain 

ds2=  {4cos2[(y-^)/2]cos2[(\(/+^)/2]}-1S(-d£2+  dy2),  and  (2.115) 

{(cos£-cosy)/(cos^+cosv)}  =  (2/Ac1)2H/S  .  (2.116) 

The    points    at    r=0    correspond    to    \|/  =  -Jt/2  or  \j/  =  7t/2;    those    at    r=°° 
correspond   to  cos^=-cos(y)tanh(7t/2/3)  .  The  Penrose  diagram  for  this 
space-time  is  shown  in  figure  9,  page  71. 

In  the  exceptional  case  y=A=0,    H=0  and  thus 

ds2  =  2drd(3  .  (2.117) 

All  null  geodesies  in  this  metric  either  have  r=constant  and  (3  varying 
linearly  with  the  affine  parameter,  or  (inconstant  and  r  varying 
linearly  with  the  affine  parameter.  Thus  on  all  incomplete  null 
geodesies  r  approaches  zero  as  the  affine  parameter  tends  to  its 
limiting  value.   Introduce  the  coordinates   £  and  \\f  by 

r  =  tan{(\|H-5)/2}  ,  (2.118) 

P  =  tan{(V-$)/2}   .  (2.119) 

Substituting   into   equation   (2.117)   results   in 
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ds2=  {2cos2[(\|/+£)/2]cos2[(y-£)/2]  yl{-d^2  +  d\j/2}  .  (2.120) 

The  points    at  \\r+l,=0  correspond  to  r=0;  those  at  \|/+£=7i   correspond  to 
r=°°.  The  Penrose  diagram  for  this  space-time  is  in  figure  10,  page  72. 

In  the  last  three  exceptional  cases,  H  is  nowhere  zero  and  thus 
the  (P,  r)  coordinate  patch  covers  the  entire  space-time.  In  the  cases 
where  H>0  introduce  the  coordinates  £,  and  \j/  by 

P  =  tan{(v-0/2}  ,  (2.121) 

2/H"1dr=  sin(v){cos[(v+|)/2]cos[(v-^)/2]}_1.  (2.122) 

In  cases     where  H<0,  introduce  the  coordinates  \|i  and  £  by 

P  =  tan{(^-V)/2}  ,  (2.123) 

2/H_1dr=  sin(v){cos[(^+\|/)/2]cos[(^-v)/2]  }_1.  (2.124) 

In  both  cases  the  metric  (2.101)  takes  the  form 

ds2=  {4cos2[(^+\|/)/2]cos2[(^-\j/)/2]  }_1IHI  {-d^2+dy2}  .       (2.125) 

In  the  exceptional  case  y=0,  A>0,  H=  -At^/3<0.  Using  equation 
(2.124)  ,  obtain 
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(6/Ar)  =  sin(\j/){cos[(^+v)/2]cos[(^-v)/2]}-1.  (2.126) 

Solving  for  r  and  substituting  into  (2.125)  results  in 


ds2  =  (3/A)  {-d%2  +  d\(/2}  /sin2^  .  (2.127) 


The  points  at  £,+\\f=K  or  £-\j/=ti  correspond  to  r=0;  those  at  ^=0   correspond 
to  r=°°.  The  Penrose  diagram  for  this  space-time  is  figure  11,  page  73. 

In  the  exceptional  case  y>0,  A=0,  H=y/r>0.  Using  equation   (2.122) 
results    in 

v2/y=    sin(V){cos[(VH)/2]cos[(V-^)/2]}-1.  (2.128) 

Solving  for  r  and  subtituting     into  (2.125),  obtain  the  metric 

ds2={16Y"1sin(\l/)cos3[(\(/+^)/2]cos3[(v-^)/2]  }-1/2{-d^2  +  d\j/2}.  (2.129) 

The  points  at  y=0  correspond  to  r=0;  those  at  y+£,=K  or  \|/-£=7i  correspond 
to  r=°°.  The  Penrose  diagram  for  this  space-time  is  figure  12,  page  74. 

Finally,  in  the  exceptional  case  y<0,  A=0,  H=y/r<0.   From  equation 
(2.124)  ,  arrive  at 

r2/y  =  sin(v){cos[(4+v)/2]cos[^-\|/)/2]}-1.  (2.130) 
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As  before,  use  this  and  equation  (2.125)  to  get 

ds2=  {16Y"1sin(^)cos3[(\(/+^)/2]cos3[(v-^)/2]  }-1/2{-d£2  +  d\j/2}  .(2.13  1) 

The  points  at  £=0  correspond  to  r=0;  those  at  £+\|/=-7c  or  t,-\\r=-K 
correspond  to  r=°° .  The  Penrose  diagram  for  this  space-time  is  figure 
13,  page  75. 

There  are  many  possibilities  for  the  global  structure  of  the 
planar  wall  space-time.  The  Penrose  diagram  for  each  of  the  manifolds 
M+  and  M"  consists  of  a  piece  of  one  of  the  Penrose  diagrams  of  figures 
8-13  bounded  by  a  timelike  curve.  The  Penrose  diagram  of  the  wall 
space-time  consists  of  the  Penrose  diagrams  of  M+  and  M"  with  the 
boundaries  identified  at  corresponding  values  of  r.  Figure  14,  page  76, 
is   the  Penrose  diagram   for  one  possible   wall   space-time.   In   this   case 

A+>0,  y+<0>  A_>  0,  and  y_>  0.  The  paths  of  the  wall  in  M+  and  M"  are 

identified. 


70 
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Figure  8.  The  Penrose  diagram  for  the  plane-symmetric 
space-time  in  the  case  y<0,  A>0. 
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Figure  9.  The  Penrose  diagram  for  the  plane-symmetric 
space-time  with  y>  0,  A  >  0.  The  inner  diagonals  correspond 
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Figure   10.  The  Penrose  diagram  for  the  planar 
symmetric   space-time  with    y=A=0. 
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Figure    11.   The   Penrose   diagram   for  the   plane-symmetric 
space-time  with    "pO,  A=0. 
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Figure    12.   The   Penrose   diagram   for  the   plane-symmetric 
space-time  in  the  case  y  >  0,  A  >  0. 
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r  =  0 


Figure    13.   The  Penrose  diagram   for  the  plane-symmetric 
space-time  in  the  case    y<0,  A=0. 
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Figure   14.  The  Penrose  diagram  for  a  possible  planar 
symmetric  wall  space-time.  A+>0,  y+<  0,  A_>  0,  y_>  0. 

The  dotted   lines   are  identified. 


CHAPTER  3 

THIN  WALLS  AND  PLANE-SYMMETRIC  FLUIDS 

WITH  PRESSURE  EQUAL  TO  ENERGY  DENSITY 


Introduction 

In  this  chapter  the  motion  of  a  thin  wall  is  treated  in  general 
relativity  for  the  case  where  the  regions  on  either  side  of  the  wall  are 
filled  with  fluid.  The  fluid  is  taken  to  be  irrotational  and 
plane-symmetric,  with  the  pressure  p  equal  to  the  energy  density  p. 
As  before,  the  wall  is  assumed  to  have  the  equation  of  state  x  =  To  , 
where  0  <_  r  <_  1  ,  x  being  the  tension  in  the  wall  and  o~  its 
energy-density.  Unlike  the  treatment  in  chapter  two,  the  wall 
space-times  are  taken  to  be  completely  reflection  symmetric  (as  in 
references  5,6,  and  7)  in  order  to  make  the  equations  tractable. 
Mathematically,  this  means   K  l+=-K  „w.      A  class  of  static    space-times 

is  found  to  permit  oscillatory  solutions  whenever  T>l/4.  In  a 
separate  class,  which  includes  a  Robertson-Walker  space-time,  the 
walls  simply  translate  through  space  without  turning  points.  The 
treatment  extends  to  fluid  space-times  previous  results  involving  thin 
walls  in  vacuum,  as  studied  in  references  5  and  6.  The  global  structure 
of  the  wall  space-time  is     discussed  in  the  final  section  of  the  chapter. 
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Taub-Tabensky  Space-Times 
To  derive  the  metric,  start  from  Einstein's  equation,  which  is 

Rab  -  (1/2)  Rgab    =    KTab      with  (3.1) 

Tab  =  (p  +  p)uaub  +pgab  ^  (32) 

where      p,   p    ,    and    u      are     the     energy     density,     pressure,      and 
four-velocity    of   the    ambient   fluid,    respectively.    Outside    the    wall    the 
scalar  field    <|>    makes  no  contribution  to  Tab    since    <(>  is    approximately 

constant  and  V((|))  is  taken  to  be  vanishingly  small. 

Tabensky  and  Taub  solved  equation  (3.1)  under  the 
assumption  that  the  energy  density  p  equaled  the  pressure  p,  and  that 
the  fluid  was  irrotational.  With  these  conditions  it  follows  from  the 
conservation   of  Tab  that  there  exists  a  function  F  such  that 

ua=   VaF/(-VcFVcF)(1/2>  ,  (3.3) 

and   that  equation   (3.1)   becomes 

Rab  =  2KPuaub  =  VaFVbF.  (3.4) 
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The  pressure  and  energy  density  of  the  fluid  is  given  by 

2Kp  =  2Kp  =  -VCF  VCF  ,  (3.5) 

while  the  equation  of  motion  is 


VcVcF  =  0.  (3.6) 


With  the  above  equations  and  the  assumption  of  plane 
symmetry,  one  can  proceed  as  in  Tabensky  and  Taub,  and  obtain  the 
solution 

ds2  =  +  f1/2exp(Q)(  -dt2  +  dz2  )  +    t(  dx2  +  dy2)  ,         (3.7) 

where   Q   is  an  arbitrary  line  integral  in  the  t-z  plane, 

Q  =  /t{(  F,t2  +  F,z2)dt  +  2F,tF,zdz}  ,  (3.8) 

with    comma    designating    partial    derivative.    In    these    coordinates    F 
satisfies 


F'zz  =  F<tt  +  t~lF'r  (3-9> 
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Note  that  either  t  or  z  can  act  as  a  time  coordinate,  depending 
on  the  choice  of  sign  in  equation  (3.7).  The  choice  of  solution  in  (3.9) 
determines    the    metric. 

Derivation  of  the  Jump  Conditions 

The   Tabensky-Taub   space-time   metric    takes   the   form 


ds2  =   eH(t,z)(-dt2  +  dz2)  +  t(dx2  +  dy2)  ,  (3.10) 


where    e=  +  1    and  H  =  O  '   ^exp(Q),  with  Q  as  defined  in  the  previous 

section.   Specify   the   hypersurface   Z  by  z=P(t).   Then  in  the  coordinates 
(  t,  z,  x,  y)    a  normal  covector  to  E  is  given  by 

Na  =  (-P,t,    1,0,0).  (3.11) 

Thus   the   unit  covector  na  and  the  T-fluid   four-velocity      v     are  given 
by 

na   =  jQ-(1/2)(  -P,t  ,    1,0,0),      and  (3.12) 

va  =  8Q-(1/2)(    1  ,  -P,t,0,0    )  ,  (3.13) 
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where  Q  =  eH_1[-(P,t)     +  1]  ,  the  norm  of  Na  ,  and  5    and  j  are  ±1  , 

corresponding   to  the  two  possible  directions   of  these  covectors.   Using 
(3.12)  and  (3.13)  ,  one  can  write 


Sab=  -vavb+   nanb     +   ^ab  •  (3-14> 

just  as  in  equation  (2.11).  Here,  qaj)dxadx'3  =  t(  dx2  +  dy2)  .  The 
intrinsic  metric  of  the  wall  will  then  be  hab  =  -vav^  +  qab  .     With  this 

expression  the  extrinsic  curvature  can  again  be  computed  from 
equation    (2.13).    A    straightforward    computation    yields 

£n(va)=  JEH-1{(Q-1/2P,t),t  +  (Q-1/2),z)va,  (3.15) 

£n(qab)  =  )cH-1r1(Q-(1'2)pft)qab  ,  (3.16) 

and    hence 

Kab  =  -vavbeH_1(  L't  +  J(Q"1/2)>Z)  +  d^^eH^t^L  ,    (3.17) 

where  L  =  jQ"1/2P,t  .  Equations  (1.10),  (3.14),  and  (3.17)  give  two 
equations  for  Kab  +  -  Kab  ,  one  in  terms  of  the  stress-energy  of  the 
wall,    and    the    other    in    terms    of    the    purely    geometric    quantities    of 
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equation    (3.17).    Equating    the    parts    tangential    and    orthogonal    to    va 
results  in  the  jump  conditions: 


(l/2)eH+'1f1L+  -  (l^eH/V1!^  -(1/2)ko  ,  (3.18) 


-eH+"1(  L+,t  +  j+(Q"1/2),z)  +  eUA(  L.,t  +  j.(Q/1/2),z)  = 

=  (r-l/2)Ko.  (3.19) 

From    vava    =   -1    and   the   above   definitions    a   useful   expression 
for  L  can  be  found: 

L  =  sH((vaVat  )2    -eH_1)1/2-  (3.20) 

Here  s=±_l,  and  will  be  interpreted  in  the  following  sections. 
Substituting   into  equation   (2.10)  results   in 

(l/2)er1s+((vaVat)2  -  eH+"1)1/2  -(l/2)ef1s.((vaVat)2  -  eH/1)1/2= 
=  -(1/2)kcj    .  (3.21) 

In  the  next  two  sections,  the  equations  developed  thus  far  will  be  used 
to  study  the  motion  of  thin  walls  in  two  particular  Tabensky-Taub 
space-times. 
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Static  Space-Time 

Choose  F=  az  ,  a  a  constant,  for  the  solution  of  equation  (3.9), 
and  let  e  =-1  in  the  metric  (3.10),  so  that  the  four-velocity  of  the 
ambient  fluid  will  be  timelike  and  at  rest  with  respect  to  these 
coordinates.  Using  equations  (3.7)  and  (3.8)  results  in  the  metric: 

ds2  =  f  1/2exp(a2t2/2)(-dz2  +  dt2)  +  t(  dx2  +  dy2)  .  (3.22) 

Note  that  z  is  the  time  coordinate  here,  and  t  a  space 
coordinate,  so  the  metric  is  static.  The  energy  density  and  pressure  of 
the  fluid  is 

p  =  p  =  (l/2)K_1a2t1/2exp(-a2t2)  ,  (3.23) 

which    is    smooth    for    all    t   >    0.    Before    invoking    the    assumption    of 
reflection   symmetry,   it  will   be  necessary   to   interpret  the   quantity      s 
in  (3.20),  and  then  examine  (3.21)  to  determine  what  values  of  s+  and 

s_  might  be  allowed.  For  e=-l  it  follows  that 


naVat  =  n°=  -sILIH"1,  (3.24) 


which    says    that,   in    the   direction    of   na,  t    decreases    when    s=l,    and 
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increases  when  s=-l.  Looking  then  at  equation  (3.21),  it  is  clear  that  if 
the  energy  density  of  the  wall,     a  ,  is  to  be  positive,     either  s+  =  l    and 

s_=-l,    or  s+  =  s_=l    with   the  additional  proviso  that  H+      >  H        .  The 

assumption  of  reflection  symmetry  leaves  only  the  former  possibility. 
Hence  in  a  direction  perpendicular  to  the  wall  and  pointing  into  M+  ,  t 
is  decreasing,  and  area  elements  are,  from  (3.22),  shrinking.  M+,  as  in 
chapter  two,  is  designated  an  "interior."  Note  that  na  points  from  M"  to 
M  +  ,  so  the  same  calculation  in  M"  must  be  done  using  -na  .  This 
means    that    s _=-l     also    corresponds    to    an    interior.    Thus    the    wall 

space-time  will  consist  of  two  interiors  glued  together  at  the  wall. 

In    view    of    the    above    considerations,    it    is    now    possible    to 
specialize  equations  (3.18)  and  (3.19)  to  the  present  case,  Ka^+=-Ka^  : 

-H+_1  t_1L+  =  -ko/2  ,  (3.25) 

2H+"1L+,t  =  (2r  -  1)ko72  .  (3.26) 

Equations  (8.4)  and  (8.5)  can  be  integrated,  yielding 

L+  =  C0t(2rM)/2,  (3.27) 

and  so  the  energy  density  of  the  T-fluid,  from  (3.18),  is 


8  5 


a  =  2K_1C0tr_1exp(-a2t2/2)  .  (3.28) 


From  (3.18)  ,  (3.20)  ,  and  (3.23)  comes  the  equation  of  motion 

{  (vaVat)2  +    H"1  }1/2    =  C0trexp(-a2t2/2)  .  (3.29) 

Squaring    both    sides    of   this    equation    leads    to    the    pseudo-Newtonian 
energy    equation 


(vaVat)2   +    V(t)    =    0  ,  (3.30) 


where 


V(t)  =  t1/2exp(-cc2t2/2){  1  -  C02t(4rM)/2exp(-a2t2/2)  }  .    (3.31) 

The  effective  energy   of  the  T-wall  is   zero   in   these  coordinates.    The 
motion  of  the  wall  can  now  be  completely  analyzed  for  every  V. 

Case  1:  r<  (1/4) 

Write  (3.31)  as  V(t)=t2rexp(-a2t2){t(1_4r)/2-  C02exp(-a2t2/2)}. 
Note   that  the  first  multiplicative  factor  is   positive  definite  for  t  >  0, 
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hence  acts  only  as  an  amplitude.  Focus  the  attention,  therefore,  on  the 
second  term,  {t(1_4r)/2-  C02exp(-a2t2/2)}.  For  small  enough  t,  this 
term  is   strictly  negative,  as  long  as  CQ   is  non-zero  (as  it  must  be,  by 

(3.23)).   As   t  increases,   t'  /'     increases    monotonically,    while    the 

exponential  term  similarly  decreases.  Therefore  there  will  be  a  point 
tmax     at  which  V(t)  ceases  to  be  negative  and  becomes  positive  for  all 

t  >  tmax    .  The  possible  motion  of  the  wall  is  simple:  it  emerges  from 

t=0,  expands  to  tmax   ,  and  returns. 

Case  2:    r=  (1/4) 

In  this  case,  (3.31)  becomes 

V(t)  =  t1/2exp(-a2t2){  1  -  C02exp(-a2t2/2)}.         (3.32) 

If  ICQI  <  1,  no  motion  is  possible.  If  ICQI  >  1  ,  then  V  will  be  negative 
on  (0,  tmax)  for  some  tmax  ,  and  positive  on  (tmax  ,  "),  as  before.  The 
wall  expands  to  tmax   ,  and  then  contracts. 
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Case  3:  (1/4)  <  r 

This    case    includes    r=l,    corresponding    to    domain    walls.    Set 

B(t)={l-C02t^4r"1^2exp(-a2t2/2)},  the  second  multiplicative  factor  on 

the  right-hand  side  of  equation  (3.31).  Note  that  B(t)  has  one  critical 
point    for    t    >    0,    which    clearly    corresponds    to    a    minimum.    It    then 

follows  that  if  the  expression   {[(4r-l)/4]"1ln(C02)  +  ln[(4r-l)/2a2]}  is 

less  than  one,  no  motion  is  possible,  while  if  it  is  greater  than  or  equal 
to  one  there  is  a  point  of  stability  around  which  the  wall  may  oscillate, 
or  remain   at  rest.   Hence  in   this   case   the   motion   of  the   wall   depends 

critically  on  the  choice  of  the  parameters  CQ,  a,  and  r. 

Dynamic  Space-Time 

Choose   F=b   ln(t),      b      a   constant,   for  the   solution   of  equation 
(3.9),  and  let    e=  +1    in  the  metric  (3.10).  This  gives 

ds2  =  tm(-dt2+  dz2)  +  t(dx2  +  dy2)  ,  (3.33) 


where  co  =  b    -1/2.  The  energy  density  and  pressure  of  the  fluid  are 


p  =  p  =  (l/2)K"1b2t"(£0+2)   .  (3.34) 
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Again,  the  ambient  fluid  is  at  rest  relative  these  coordinates. 
The  metric  is  singular  at  t  =  0,  and  the  wall  space-time  consists  of  two 
interiors  glued  together,  just  as  in  the  previous  case.   Proceeding  in  the 

same   way   as   before   yields   L+  =  -C    t       '    .  The    energy  density    of  the 
T-fluid  is     given  by 

a=2K-1Cotr-(0-3/2.  (3.35) 

An   explicit   expression   for   the   ordinary    velocity,    dz/dt,   can    be   easily 
found: 


dz/dt  =  +C0tP(  C02t2P  +  1  )_1/2  ,  (3.36) 


where  p  =  T-(b  /2)  -  (1/4)  .  While  this  expression  can  be  exactly 
integrated  only  for  special  values  of  the  parameters,  it  is  nonetheless 
clear  by  inspection  what  the  possible  motions  of  the  wall  must  be. 

Case  1:    p  >  0 


Here  the  wall  starts  from  rest  at  some  zQ   and  accelerates   to  the 
speed  of  light  in  the  limit  as      t      goes  to  infinity. 
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Case  2.   p  =  0 

In  this  case  the  wall  proceeds  at  some  uniform  velocity  in  either 
the   positive   or  negative   z-direction. 

Case  3:  (3  <  0 

The  wall  moves  initially  at  the  speed  of  light  and  slows  to  zero 
in  the  limit  as    t    goes  to  infinity. 

Case  4.       b=(3/2)1/2 

This  corresponds  to  a  Robertson-Walker  metric,  as  can  be  seen 
by  setting  T=(2/3)r'  ,  which  puts  the  metric  in  the  usual  form, 
ds2=-dT2+a(T)(dx2+dy2+  dz2  ).  Here  (3  =  T-l,  which  means  that  domain 
walls  will  translate  at  uniform  velocity  through  the  space-time.  As  the 
wall  expands,  space  is  opened  up  and  fluid  is  created  at  the  expense 
of  the  energy  density  of  the  wall. 

Global  Structure 

As  discussed  previously,  for  the  space-times  considered  here 
one  must  glue  together  two  interiors.  In  both  cases,  there  is  a  real 
singularity   at   t=0,    time-like   for   the    static    space-times    and    space-like 
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for   the   dynamic    space-times.   These   space-times   are   topologically   R   . 
To  draw  the  Penrose  diagrams,     consider  the  two  dimensional  metric 


ds2=  eH(t,z)[-dt2+dz2]  .  (3.37) 


Each  point  of  the  diagram  will  correspond  to  a  2-plane.  The  null 
geodesies  issue  from  the  singularity  at  t=0  and  extend  to  arbitrarily 
large  parameter  values.  Since  t  is  strictly  positive  in  these  space-times, 
the  Penrose  diagrams  will  be  very  similar  to  the  diagram  for 
Minkowski  space  in  spherical  coordinates  ,  with  the  exception  that 
t=0  will  correspond  to  a  real  singularity,  rather  than  simply  a 
coordinate    singularity.    The   desired    coordinate   transformation    is 

t=  tan((t*+z*)/2)  +  tan((t*-z*)/2)  ,  (3.38) 

z  =  tan((t*+z*)/2)  -  tan((t*-z*)/2)  .  (3.39) 

This    transformation    yields    the    two-metric 

ds2=  eHsec2((t*+z*)/2)sec2((t*-z*)/2){-dt*2+dz*2}  .  (3.40) 

Figures  15  and  16,  pages  92  and  93,  are  Penrose  diagrams  for  the 
static  and  dynamic  cases,  respectively,  each  with  a  sample  trajectory. 
To  construct  the  Penrose  diagram  for  the  T-wall  space-times  studied 
here,  cut  along  the  trajectory  of  the  wall,  taking  the  piece  in  which  the 
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area  elements  are  shrinking  in  the  direction  normal  to  the  wall.  This 
piece  corresponds  to  the  interior.  Gluing  this  piece  to  its  mirror  image 
gives   the   desired   diagram. 


92 


t=0 


t=oo 


Figure  15.  The  Penrose  diagram  for  static  Tabensky-Taub 
space-time,  with  sample  trajectory. 
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t=oo 


t=oo 


z=oo 


t  =  0 


Figure  16.  The  Penrose  diagram  for  dynamic  Tabensky-Taub 
space-time,  with  sample  trajectory. 


CHAPTER  4 
THIN  WALLS  IN  DUST  UNIVERSES 

Introduction 


In  this  chapter  the  motion  of  a  thin  wall  is  treated  in  general 
relativity  for  the  case  where  the  regions  on  either  side  of  the  wall  are 
filled  with  dust.  As  in  chapter  three,  the  space-time  is  taken  to  be 
plane    and    reflection-symmetric,    Ka^  +  =  -K&^)_,   and  dependent   only   on 

time.  As  before,  the  wall  is  assumed  to  have  equation  of  state  x  =  To  , 
0<T<1,  where  x  is  the  tension  and  o*  is  the  energy  density  of  the  wall. 
Two  classes  of  exact  solutions  to  Einstein's  equation  are  presented, 
both  corresponding  to  a  plane-symmetric  dust  with  cosmological 
constant.  To  make  the  equations  tractable,  the  arbitrary  function  of  z 
in  these  metrics  is  taken  to  be  a  constant,  so  the  z-dependence  is 
eliminated.  Using  these  metrics,  classes  of  space-times  are  found  in 
which  the  thin  wall  undergoes  translation  without  turning  points: 
initially  comoving,  the  wall  accelerates  to  some  maximum  velocity,  and 
returns  to  comoving  in  the  limit  as  time  goes  to  infinity.  The  wall 
tends  to  a  limiting  value  of  the  space  coordinate.  Other  classes  are 
found  in  which  the  wall  exits  from  the  initial  singularity  at  high 
velocity,  and  again  returns  to  comoving  with  the  ambient  fluid. 
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Plane-Symmetric  Dust  Space-Times 
Einstein's  equation  is  given  by 

Rab-   (l/2)Rgab  +  Agab  =  KTab  .  (4.1) 

Here     T   b  =  puauu  ,  p    the    energy    density    of    the    dust    and    ua    its 

four-velocity.  Using  co-moving  coordinates  and  starting  from  the 
metric 

ds2=  -e2v(t'z^dt2  +  e2^1'2^2  +  Y(t)2[dx2  +   dy2]  (4.2) 

it  is  possible  to  derive  an  exact  solution  of  (4.1),  which  bifurcates 
when  one  of  the  parameters  goes  to  zero.  V  Ta  =  0  gives  v,z=0.  The 
field  equations   are     : 

G0°  =  -2Y_1e"2v{  Y,tX,t  +  2~1(Y,t)2Y-1 }  +  A  =  -kP  ,  (4.3a) 

Gj1  =  -2Y"V2v{  Y,tt-Y,tv,t  +  2"1(Y,t)2Y"1}  +A  =  0,       (4.3b) 
G22=G33=-e-2v{?i,tt  +(?,,t)2^,tv,t  +Y-!Y,tt  +  Y^Y^X,,  -  v,t)}  + 
+  A  =0  .  (4.3c) 
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Take  Y(t)  as  the  time  coordinate,  i.e.  Y=t  .  Then  equation  (4.3b)  can  be 
integrated    by 

e"2v  =    (A/3)t2  +  (c/t)  ,  (4.4) 

where  c  is  an  arbitrary  constant  of  integration.  Next,  following  the 
general  method  of  Kramer  et  al.  ,  set  e  -e"vV(t,z).  Inserting  this  and 
equation  (4.4),  with  c*0,  into  equation  (4.3c)  yields 

V(t,z)=  A(z)tev  +  B(z)   ,  (4.5a) 

where    A(z)    and    B(z)    are    arbitrary    functions    of    integration,    one  of 

which    can    be    eliminated    by    a    transformation    of    the    coordinate  z. 

Taking  c=0  in  equation  (4.4)  and  proceeding  as  before  gives  an 
alternate    solution, 

V(t,z)  =  A(z)e3v  +  B(z)  .  (4.5b) 

Note  that  taking  the  limit  as  c-»  0  in  (4.5a)  results  in  V  becoming  a 
function  of  the  coordinate  z  alone,  whereas  in  (4.5b)  V  is  a  function  of 
both  t  and  z.  In  the  limit  as  A->0  ,  the  space-times  corresponding  to 
(4.5a)  go  over  to  the  usual  plane-symmetric  dust  solutions,  whereas 
the  space-times  corresponding  to  (4.5b)  link  to  plane-symmetric 
vacuum    solutions.    This    can    be    best    appreciated    by    obtaining    the 
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energy  density  p   from  equation  (4.3a),  which  is: 

kP  =   (2/t)[  (A/3)t2  +  (c/t)  ]  V^V1.  (4.6) 

This  equation  shows  that,  in  the  c=0  space-times,  the  energy  density  of 
the  dust  is  coupled  directly  to  the  vacuum  energy,  in  contrast  to  the 
c*0  case.  The  metrics  in  the  two  cases  are  explicitly: 

ds2=-[(A/3)t2+  c/tl^dt2  + 
+  [(A/3)t2  +  c/t]{At[(A/3)t2+c/tT1/2+B}2dz2  + 

+  t2[dx2+dy2]  ,  and  (4.7a) 

ds2=-(3/A)f  2dt2  +  (A/3)t2[(A/3)"3/2At"3  +  B]2dz2+ 

+  t2[dx2+dy2]  .  (4.7b) 

The  metric  coefficients  in  (4.7a)  resemble  those  of  the 
plane-symmetric  Schwarzschild-De  Sitter  metric,  with  t  and  z 
reversing  roles.  These  solutions,  apparently  new,  are  interesting, 
particularly  in  that  they  both  approach  Robertson-Walker  space-time 
in  the  limit  of  large  t,  when  B  is  a  constant.  They  may  be  taken  as 
non-homogeneous  dust  universes  which,  for  appropriate  choice  of  A 
and  B,     approach  a  standard  model  as  t  becomes  large. 
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Derivation  of  the  Jump  Conditions 

Next,  we  turn  to  the  equations  of  Israel  giving  the  jump  in  the 
extrinsic  curvature  across  a  thin  shell  of  stress-energy.  The  calculation 
proceeds  precisely  as  in  previous  chapters.   Let 


ds2  =  -f(t)dt2  +  h(t,z)dz2  +  t2[dx2  +  dy2]  .  (4.8) 


Specify  the  hypersurface  £  by  z=P(t),  as  in  chapter  three.  In 
the  coordinates  (t  ,z  ,x,  y)  a  normal  covector  to  £  is  given  by 
N=(-P,t  ,1   ,0,0),  where  the  comma  denotes  partial  derivative.  Then 

na    and    va    are  given  by 

na  =  jQ"1/2(  -P,t,  1  ,  0  ,0  )  ,     and  (4.9) 

va  =  5Q"1/2(  fl/2h-l/2  ^      _rl/2hl/2p>t     ^     0,0),         (4.10) 

1  1  9 

where     Q  =[h      -  F   (P,t)   ]  is  the  squared  norm  of  N&  ,  and  j  and  5  are 

±1,    corresponding    to    the    two   possible   choices    of   direction    for    these 
covectors.   Equation   (2.11),      g  ^  =  -vaVj?  +    n&n^  +  qa^  ,  holds  just  as 

before,   with   qa^dxadxb  =  t  (dx     +  dy   )  .     The  intrinsic  metric  of  the 

wall  is  again  ha^  =  -vav^  +  q  l    ,   and   the  extrinsic   curvature  can   be 
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computed     from        equation     (2.13),        Kab     =    (l/2)£n  (hab).     This 
computation    results    in 

Kab=  -(fh)"1/2!  L,t  +  (P,tr1L>z}vavb   +     {(fhr^t^LJq^  ,    (4.11) 

where     L=  j(P,  t)(f  ^Q"1)1/2.   Inserting  (4.11)  into  the  left  hand   side 
of  (1.10)  and  equating  parts  parallel  and  normal  to     va   gives  the  jump 

conditions,     which,     specialized     to     the     assumptions     of    reflection 
symmetry   and   dependence   only   on   time,   yield: 

(f+h+)"1/2L+    =  -ica(2r-l)/4  ,        and  (4.12) 

(f+h+)"1/2t_1L+    =-ko/4  .  (4.13) 

Eliminating     a    and     integrating  these  equations  yields 


L+  =  -C0t2lM.  (4.14) 


Now      (vaVa  t)2  =    f1    +  (fh)-1L  +  2    ,  and  (vaVaz)2  =  h-2L  +  2  ,    so 

together  with  (4.14)  it  is  possible  to  obtain  an  explicit  expression  for 
the  ordinary  velocity  of  the  wall  as  a  function  of  time,  just  as  was  done 
in  Chapter  Three: 
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Idz/dtl     =    C0f1/2h"1/2[ht2-4r  +    C02]"1/2.  (4.15) 

Wall  Dynamics 

Using  equation  (4.15),  it  is  now  possible  to  analyze  the  motion 
of  the  wall.  Several  cases  will  be  considered,  with  A,  B,  c,  and  A  taken 
to   be   non-negative   constants. 

Case  l:c>0.  A>0.  B>0.  A*0 

Here,   equation   (4.8)   becomes 
ds2  =  -fdt2  +  f[  Atf1/2  +  B]2dz2  +  t2[  dx2  +  dy2]  ,  (4.16) 

where  f=[  (A/3)t2  +  c/t  ]_1,  while 

kP  =  3Acf2{At  +  B  [  (A/3)t2  +  c/t  ]1/2}_1.  (4.17) 

As  t-»  0  ,  the  energy  density  Kp->°°,  while  f~t  and  h~t  ,  so  that 
ldz/dtl~t2r+1'2,  from  (4.15),  which  goes  to  zero  in  the  limit  as  t-»0.  As 
t->°°,  Kp->0,  f~t~2  and  h~t2,  and  (4.15)  gives  Idz/dtl-  t2r  ,  so  again  the 
velocity  goes  to  zero  in  the  limit.  Since  there  are  no  turning  points,  the 
gross  features  of  the  motion  of  the  wall  are:  the  wall  exits  from  the 
initial    singularity    while    comoving    with    the    ambient   fluid,    accelerates 
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to  some  maximum  value,  then  decelerates  back  to  the  comoving  state. 
In  these  coordinates,  it  also  follows  that  the  value  of  z  at  which  the 
wall   returns   to   comoving   is   asymptotically   approached. 

Case  2:  c>Q.  A>0  .  B>0.  A  =  0 

With    vacuum   energy   equal    to    zero,    the    metric    takes    the   form 
(4.16)  with  f=t/c  .  Explicitly,  this  is 

ds2=  -(t/c)dt2  +  (c/t){Ac"1/2t3/2+B}2dz2  +  t2(dx2  +  dy2).    (4.18) 

The  energy  density  is  given  by 

Kp  =  3Act"2{At+B(c/t)1/2}_1,  (4.19) 

which  diverges  as  t-»0  and  falls  off  to  zero  as  t->°°.  As  t->0  ,  h~t  ,  so 
clearly  (4.15)  gives  Idz/dth t  if  l-4r>0,  and  ~t2r  +  1/2  if  l-4r<0.  In 
either  case,  the  velocity  goes  to  zero  in  the  limit.  As  t-*00,  h~t  ,  and 
ldz/dtl~t  J/z,-»0  for  the  values  of  r  under  consideration.  Hence  the 
qualitative  behavior  is  the  same  as  in  case  1  as  long  as  T<3/4. 
Otherwise,  the  wall  will  not  asymptotically  approach  a  limiting  value 
of  the  z-coordinate  as  t->°°. 
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Case  3:  c>0  .  A>0.  B=0.  A*0 

This     case     corresponds     to     a     class     of     Robertson-Walker 
space-times.   The  metric   and  energy  density   are  given   by 

ds2=  -[(A/3)t2  +  c/t]_1dt2  +  t2[  dz2  +  dx2  +  dy2]  ,  and         (4.20) 

Kp  =  3ct"3  .  (4.21) 

Setting   t=ccsinh  '3(pT),  with  a  and  p  constants,  converts  (4.20)  into  the 
usual  form:  ds2=-dT2+a(T)(dz2+dx2+dy2).  In  this  case,  as  t-»0,  f~t  and 

h~t2,  so    from    (4.15)     ldz/dtl~  t"  ^2.  For  t^<»,f~t"2    and   h~t2,    so 

IV  4 
ldz/dtl~t  .  It  is  easy  to  see  that  for  all  values  of  T   between  zero  and 

one,  the  wall  exits  the  singularity  at  high  velocity  and  falls  off  to  zero 

velocity  in   the   limit  as      t-»  * ,    approaching    some    finite    value    of    z 

asymptotically. 

Case  4:  c>Q.  A>0.  B=0.  A=0 

As      in     case     three,      these     space-times      are     obviously     of 
Robertson-Walker   type.    The   metric   and   energy   density   are 


ds2  =  -(t/c)dt2  +  t2(dz2  +  dx2  +  dy2),     and  (4.22) 
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kP  =  3cr3.  (4-23) 

Here,  f~t  and  h~t2,  so  ldz/dtl~t~1/2  as  t->0.  As  t->«,  ldz/dtl~t2r~5/2.  As 
in  case  two,  the  wall  exits  at  high  velocity  from  the  initial  singularity 
and  gradually  falls  to  zero  velocity,  though  it  remains  at  bounded  z  in 
these  coordinates  only  when  r<3/4. 

Case  5:   c=0.  A*0 

In  this  case,  the  metric  (4.7b)  must  be  used,  since  c=0  in  (4.7a) 
gives  p=0.  In  the  case  of  (4.7b),  the  metric  and  energy  density  of  the 
dust  are  given  by 

ds2  =  -(3/A)t~2dt2  +  (A/3)t2{(A/3)-3/2At"3+B}2dz2  + 

+  t2{dx2  +  dy2},    and  (4.24) 

kP  =  -2AA(A  +  Bt3(A/3)3/2}_1.  (4.25) 

A  and  B   are  here  taken  to  be  constants,  their  precise  values  having  a 
direct   effect   on    the   permitted   domain    of   t.    Choosing    A=-B(A/3)3' 
mandates  t>l    to  avoid   a  dust  with  negative  energy  density      for  values 

less    than    t=l.    As    t    goes    to    1,    there    is    again    a    singularity,    with 

3         1 
ldz/dtl~t(t   -1)     .  The  wall  exits  at  high  velocity  and  goes  to  infinity  in 

the  limit  as  t  -*00,  while  its  velocity  falls  off  to  zero  like  t  . 
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Global  Structure 

As  in  chapter  three,  the  global  structure  of  the  wall  space-time 
consists  of  two  interiors  glued  together,  one  the  mirror  image  of  the 
other.  As  before,  begin  by  studying  the  behavior  of  plane-orthogonal 
null   geodesies.   Consider  the   2-dimensional   metric 

ds2=  -f(t)dt2  +  h(t)dz2.  (4.26) 

Then,  as  in  chapter  two,  there  exists  a  constant    co    such  that 

co  =  kb(3/az)b=  kj  =  h(dzAR).  (4.27) 

For  null  geodesies, 

0  =  -f(dt/d>.)2  +  h(dzAR)2=  -f(dtAR)2  +  h^co2,  (4.28) 

and  so  it  follows  that 

cod*.  =±(fh)1/2dt  .  (4.29) 

Hence  in  case  1  where  A,  B,  c,  and  A  are  positive  constants, 
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co  /  dX  =  t/0"  { At[(A/3)t2+c/t]'1/2  +  B}  dt  .  (4.30) 

For  small  t  ,  coA.  =  at  '  +  Bt,  a  a  constant,  plus  an  arbitrary  constant 
which  can  be  set  equal  to  zero.  So  as  X  goes  to  zero,  all 
plane-orthogonal  null  geodesies  encounter  the  singularity  at  t=0  and 
are  past  incomplete.  On  the  other  hand,  for  large  t,  (aX  =  at  +  d,  a  and  d 
constants,  so  that  as  t-»  °°  so  does  X,  and  the  plane-orthogonal  null 
geodesies  are  future  complete  and  inextendible.  Now  define  the 
coordinate  T  by 

Ts  /Qt  (f/h)1/^!.  (4.31) 

This  coordinate  transformation  will  be  used  in  each  of  the  following 
cases,  sometimes  taken  to  be  improper  in  its  lower  limit.  The 
integrand  is  positive  definite,  so  T(t)  is  a  strictly  increasing  function  of 
t.  Since  the  integrand  goes  like  t  for  small  t  and  t  for  large  t,  T  is 
bounded  above.   Set  Tmax=/0<*>(f/h)1'2dt  .  The  metric  is  given  by 

ds2  =  h(-dT2  +  dz2)  ,  (4.32) 

and  the  Penrose  diagram  is  displayed  in  figure   17,  page   108. 

Next  consider  the  case  2:  c>0,  A>0,  B>0,  and  A=0.  Define  T(t)  as 
before,  this  time   obtaining 
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T(t)  =  /0l  t{/cAt3/2+  Bc}_1dt  ,  (4.33) 

which    in    this    case   diverges    in    the    limit    as    t    goes    to    infinity.    The 
coordinate  T  puts  the  metric  in  the  form  (4.32).  Now  define 

T  =  tan{(t*+z*)/2)  +  tan{(t*-z*)/2}  ,  (4.34a) 

z  =    tan{(t*+z*)/2}  -  tan{(t*-z*)/2}  .  (4.34b) 

The  resulting  metric  is 

ds2  =  hsec2{(t*+z*)/2}sec2{(t*-z*)/2}{-dt*2  +  dz*2}.  (4.35) 

The  Penrose  diagram  is  given  in  figure   18,  page  109. 

In  case  3,  c>0,  A>0,  A>0,  and  B=0.  Change  the  time  coordinate 
with  T(t)  =  lima_>0/at  t_1{(A/3)t2  +  c/t}"1/2dt  .  The  Penrose  diagram 

is  identical  to  that  of  case  1,  given  in  figure  17,  page  108. 

In  case  4  ,  with  c>0,  A>0,  B=0,  and  A=0,  the  time  coordinate 
transformation  (4.31)  results  in  t=(cT/2)  ,  putting  the  metric  in  the 
form  ds-(cT/2)2{  -dT2  +  dz2}.  Here,  using  the  transformation 
(4.34a,b),  a  metric  of  the  form  (4.35)  is  again  obtained.  The  Penrose 
diagram  is  identical  to  that  of  case  2,     shown  in  figure  18,  p.  109. 

Putting    c=0,    A>0,    choose    metric    (4.7b),    which    gives    case    5. 
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Choosing  B=-(3/A)^'2A  makes  the  energy  density  positive  definite  for 
t  >  1,  and  places  the  singularity  at  t=l.  Other  choices  lead  to  similar 
space-times.  For  example,  it  is  possible  to  arrange  matters  so  that  the 
singularity  is  in  the  future  rather  than  the  past,  amounting  to  a 
reflection  in  time.  With  the  given  choices,  it  can  easily  be  found  that, 
for  all  null  geodesies  orthogonal  to  the  planes  of  symmetry, 


(oX  =  t  +(2t2)_1  +  cj  ,  (4.36) 


where  equation  (4.29)  was  used.  From  (4.36),  it  is  clear  that  the 
plane-orthogonal  null  geodesies  are  past  incomplete  and  future 
complete.    The    transformation    (4.31),    slightly    modified,    now    gives    the 

new    coordinate    T(t)=(A/3)/t°°  t(l-;r )~   dt.  Note   that  infinity   and   zero 

exchange  places  in  this  coordinate  transformation.  Using  the 
transformations  (4.34a,b),  it  is  then  possible  to  put  the  metric  in  the 
form  (4.35).  The  Penrose  diagram  is  given  in  figure   19,  page  110. 

The  Penrose  diagram  of  the  wall  space-time  may  be  obtained  as 
in  chapter  4,  figures  15  and  16,  by  cutting  along  the  trajectory  of  the 
wall,  taking  the  interior  piece,  and  attaching  it  to  a  mirror-image.  Such 
a  space-time  is  illustrated  in  figure  20,  page   111. 
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t  =  oo 


t=0 


Figure  17.  The  Penrose  diagram  for  a  plane-symmetric  dust 
space-time  in  the  case  c>0,  A>0,  B>0,  and  A>0.  The  same 
diagram  holds  for  case  3  (constants  as  in  case  1 ,  with  B=0). 
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t=oo 


t=0 


Figure  18.  The  Penrose  diagram  for  the  plane-symmetric 
dust  space-time  of  case  2,  with  c>0,  A>0,  B>0,  and  A=0.  The 
same  diagram  holds  for  case  4  (c,  A,  B,  as  above,  but  with 
B=0). 
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Figure  19.  The  Penrose  diagram  for  the  plane-symmetric 
dust  space-time  of  case  5,  with  c=0,  A>0,  B>0,  and  A>0. 
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Figure  20.  The  Penrose  diagram  for  the  plane-symmetric 
wall  space-time  of  case  2,  with  c>0,  A>0,  B=0,  and  A>0.  The 
left-hand  side  is  an  interior,  and  the  right-hand  side  its 
reflection.  The  two  pieces  are  identified  along  the  dotted 
lines,  the  trajectory  of  the  wall. 


CHAPTER  5 
CONCLUSIONS 


In  the  course  of  this  work  the  dynamics  of  thin  walls  have  been 
examined  in  space-times  with  non-zero  stress-energy.  Specifically,  the 
space-times  studied  here  included:  (i)  spherically  symmetric  with 
vacuum  energy,  (ii)  plane-symmetric  with  vacuum  energy,  (iii) 
plane-symmetric  fluids  with  pressure  equal  to  energy  density,  and 
(iv)  plane-symmetric  dust  with  vacuum  energy.  Each  of  these  will  be 
discussed   in   turn. 

The  dynamics  of  a  spherical  or  planar  T-wall  in  regions  of 
vacuum  energy  depend  on  several  different  parameters:  the  vacuum 
energy  A  and  the  "mass  parameter"  m  or  y  of  the  regions;  the  equation 
of    state    parameter  V  and    the   constant   cQ    which    are   related    to    the 

energy  density  and  tension  of  the  wall.  The  dependence  of  the  wall 
motion  on  these  parameters  is  fairly  complicated  and  there  are  many 
possibilities,  which  have  been  examined  case  by  case,  for  the  motion  of 
the  wall.  This  work  generalizes  the  work  of  previous  authors  who 
examined  wall  dynamics  for  various   special   values   of  the  parameters. 

In  all  cases  the  problem  of  the  wall's  motion  is  reduced  to  an 
equivalent  one-dimensional  mechanics  problem.  The  radius  R  of  the 
wall,  proper  time  along  the  wall,  and  the  function  V_(R)  play  the  roles 
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of  spatial  coordinate,  time,  and  potential,  respectively.  The  form  of  the 
potential   depends    on   the   values   of  the   parameters   A  +  ,  A  ,   and   either 

m+  ,  m_,  or  y+,  y_. 

In  the  case  of  spherical  walls  the  possible  values  of  the 
parameters  can  be  divided  into  a  generic  case  and  several  exceptional 
cases.  In  all  the  exceptional  cases  the  potential  has  at  most  one 
maximum  and  no  minima.  Thus  the  wall's  motion  has  at  most  one 
turning  point  which  is  a  point  where  dR/d>.  changes  sign.  In  the 
generic  case  for  sufficiently  "stiff"  equations  of  state,  i.e.  V  greater  than 
some  value,  the  potential  has  exactly  one  maximum  and  no  minima. 
Therefore  the  wall's  motion  has  at  most  one  turning  point.  A  domain 
wall  (r=l)  is  a  particular  example  of  a  wall  with  a  stiff  equation  of 
state.  In  the  generic  case  for  sufficiently  "soft"  equations  of  state,  i.e.  r 
less  than  some  value,  there  are  always  values  of  the  other  parameters 
for  which  the  potential  has  a  local  minimum.  Thus  the  wall's  motion 
can  have  two  turning  points,  which  means  it  can  undergo  oscillations 
about  a  radius  of  stable  equilibrium.  A  dust  wall  (r=0)  is  a  particular 
example  of  a  wall  with  a  soft  equation  of  state. 

In    the    case    of    planar    walls,    the    motion    exhibits    no    turning 

points  if  the  parameter  y+  +  y_     is  nonpositive.  When  y+  +  y    is  positive 

the  possible  values  of  the  parameters  can  be  divided  into  a  generic 
case  and  several  exceptional  cases.  In  all  of  the  exceptional  cases  the 
potential    has    at    most    one    maximum    and    no    minima.    Therefore    the 
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wall's  motion  has  at  most  one  turning  point.  In  the  generic  case  for 
sufficiently  stiff  equations  of  state  the  potential  has  exactly  one 
maximum  and  no  minima.  Therefore  the  wall's  motion  has  at  most  one 
turning  point.  In  the  generic  case  for  sufficiently  soft  equations  of 
state  there  are  always  values  of  the  other  parameters  for  which  the 
potential  has  a  local  minimum.  So,  as  in  the  spherical  case,  the  wall's 
motion  can  have  two  turning  points  and  there  can  be  oscillation  about 
a  radius   of  stable  equilibrium. 

In  contrast  to  the  results  for  static  planar  space-times  with 
vacuum  energy  ,  static  fluid  space-times  with  p  =  p  do  not  permit 
oscillatory  solutions  for  r  <  (1/4).  On  the  other  hand,  such  solutions  are 
permitted  for  r  >  (1/4),  whereas  in  the  vacuum  energy  case  they  are 
not  for  T  >  (8-3/  2)/8  .  This  reversal  appears  to  be  due  to  the 
exchanging  of  the  roles  of  t  and  z  in  this  class  of  space-times.  For  the 
dynamic  space-times,  uniform  velocity,  accelerated  motion,  and 
decelerated  motion  are  all  possible,  without  turning  points.  When  the 
parameters  are  chosen  appropriately,  the  wall  decelerates  from  high 
velocity  after  exiting  from  the  initial  singularity.  The  energy  density 
both  of  the  wall  and  of  the  ambient  fluid  fall  off  with  the  ensuing 
expansion,  as  would  be  expected  in  a  model  of  an  expanding  universe. 

The  dynamic  dust  universe  is  similar  to  the  dynamic  universe 
with  p=p  in  that  there  are  no  turning  points  in  the  cases  under  study. 
It  is  interesting  that  the  cosmological  constant  term  makes  no  essential 
difference  for  the  classes  of  parameter  values   chosen   here,   as   long   as 
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r<3/4.  For  domain  walls,  with  r=l,  the  difference  involves  asymptotic 
approach  to  a  limiting  value  of  z  in  the  presence  of  non-zero  vacuum 
energy,  as  compared  to  escape  to  infinity  in  the  absence  of  vacuum 
energy.  In  every  case  the  wall  moves  without  turning  points 
regardless  of  the  value  of  r.  It  is  interesting  to  conjecture  as  to 
whether  or  not  an  absence  of  turning  points  is  a  feature  of  dynamic 
space-times  in  general.  If  the  universe  were  closed,  such  turning 
points  would  have  to  exist;  in  open  universes,  they  wouldn't  be 
expected,  except  perhaps  locally.  The  Schwarzschild-de  Sitter  metric 
and  its  plane-symmetric  analogue  are  descriptive  of  essentially  local 
pieces  of  space-time,  whereas  the  fluid  and  dust  space-times  are  more 
appropriate  as  descriptions  of  the  universe  in  the  large.  Again,  it  isn't 
clear  what  inferences  can  be  made  without  a  complete  solution  to  the 
general    problem. 

As  seen  in  the  text,  there  are  also  many  possibilities  for  the 
global  structure  of  the  wall  space-time.  In  the  vacuum  energy  "generic 
cases",  i.e.  cases  where  the  vacuum  energies  and  mass  parameters  do 
not  have  special  values,  we  find  that  in  the  spherical  case  the  regions 
on  either  side  of  the  wall  have  a  global  structure  either  as  depicted  in 
figure  1  (if  9m2A>l),  page  40,  or  figure  3  (if  0<  9m2A<l),  page  42.  Here 
the  singularities  and  infinity  are  spacelike.  In  the  generic  planar 
vacuum  energy  case  the  global  structure  of  these  regions  is  the  sort 
depicted  in  either  figure  8(if  y<0,  A>0),  page  68,  or  figure  9  (if  y  >0, 
A>0),  page  69.  Here  again  infinity  is  spacelike,  and  the  singularities  are 
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spacelike  for  y  <0  and  timelike  for  y  >0.  These  global  structures  are 
quite  different  from  those  of  de  Sitter,  Minkowski,  and  Schwarzschild 
space-times  and  their  planar  analogues.  Thus  a  generic  T-wall  in 
regions  of  vacuum  energy  has  a  structure  very  different  from  that  of 
the  simplest  false  vacuum  bubbles  or  that  of  a  T-wall  in  vacuum.  This 
is  also  true,  to  a  lesser  extent,  for  walls  in  the  fluid  and  dust 
space-times  studied.  These  space-times,  mainly  because  of  the 
assumption  of  reflection  symmetry,  have  a  simpler  structure  than 
most  of  their  counterparts  in  the  vacuum  energy  case.  The  assumption 
of  reflection  symmetry,  while  physically  reasonable,  greatly  reduces 
the  generality.  In  chapter  three,  the  static  Tabensky-Taub  space-times 
exhibit  a  timelike  singularity,  while  the  dynamic  space-times  have  a 
spacelike  singularity.  The  dust  metrics  of  chapter  four  are  of  three 
types.  Figure  17,  page  106,  exhibits  a  spacelike  singularity  and 
spacelike  infinity,  which  is  the  general  case  with  the  parameters  c,  A, 
B,  and  A  all  positive.  Setting  B=0  doesn't  change  the  diagram.  Taking  c, 
A,  and  B  positive,  and  A=0  gives  the  diagram  of  figure  18,  page  107,  in 
which  the  singularity  is  again  spacelike.  As  before,  setting  B=0  does 
not  affect  the  diagram.  The  final  type,  with  c=0,  and  A,  B,  and  A 
positive,  corresponds  to  the  alternate  dust  metric,  (4.7b).  The  diagram 
on  page  108  exhibits  spacelike  infinity  while  the  singularity  is  null  in 
these  coordinates  .  The  corresponding  wall  space-times  are  created 
with  mirror  images  of  space-time  interiors,  for  both  the  dust  and  fluid 
cases. 
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There  remain  a  number  of  interesting  unsolved  problems 
concerning  thin  walls  in  general  relativity.  For  example,  the  dynamics 
of  a  radiating  wall  poses  the  difficult  problem  that  the  jump  conditions 
depend  on  both  a  time  and  a  space  coordinate,  so  the  equations  are 
partial  differential  equations,  not  the  ordinary  differential  equations 
treated  in  this  and  previous  work.  There  is  also  the  question  as  to 
whether  or  not  thicker  walls  would  be  expected  to  exhibit  the  same 
dynamics  as  the  thin  walls  of  Israel's  approximation,  which  are  taken 
to  have  negligible  thickness.  Another  question  concerns  the  fact  that, 
in  fluid  and  dust  universes,  ambient  stress-energy  is  created  out  of 
the  energy  density  of  the  wall  as  it  expands,  or  conversely  is  absorbed 
as  the  wall  contracts.  This  is  puzzling  and  somewhat  artificial,  since 
there  is  no  mechanism  for  it.  The  interesting  and  important  point, 
however,  is  that  Israel's  method  can  be  successfully  applied  to 
space-times  containing  stress-energy,  and  that  exact  answers  can  be 
obtained  concerning   the  dynamics  of  the   walls. 
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